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Abstract: We consider the recently developed diagrammatic approach to soft-gluon ex- 
ponentiation in multiparton scattering amphtudes, where the exponent is written as a sum 
of webs ~ closed sets of diagrams whose colour and kinematic parts are entangled via mix- 
ing matrices. A complementary approach to exponentiation is based on the multiplicative 
renormalizability of intersecting Wilson lines, and their subsequent finite anomalous di- 
mension. Relating this framework to that of webs, we derive renormalization constraints 
expressing all multiple poles of any given web in terms of lower-order webs. We examine 
these constraints explicitly up to four loops, and find that they are realised through the 
action of the web mixing matrices in conjunction with the fact that multiple pole terms 
in each diagram reduce to sums of products of lower- loop integrals. Relevant singularities 
of multi-eikonal amplitudes up to three loops are calculated in dimensional regularization 
using an exponential infrared regulator. Finally, we formulate a new conjecture for web 
mixing matrices, involving a weighted sum over column entries. Our results form an im- 
portant step in understanding non-Abelian exponentiation in multiparton amplitudes, and 
pave the way for higher-loop computations of the soft anomalous dimension. 
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1 Introduction 

The infrared singularities of scattering amplitudes [1-40] are important both for pragmatic 
collider phenomenology applications and for theoretical reasons. The singularities govern 
the structure of large logarithmic corrections to cross sections to all orders in perturbation 
theory, allowing the possibility to resum such terms, see e.g. [41-71] . They have also been 
explored in a variety of more formal contexts, such as supersymmetric gauge theories and 
quantum gravity [72-81], in order to examine the higher loop structure of those theories 
and to elucidate the relationships between them. 

It is well known that infrared (IR) singularities can be traced to the emission of soft 
and collinear gluons, where the latter contribute only in the case of massless partons. 
Consequently an amplitude may be factorised into hard, soft and jet functions, in such 
a way that the hard function, describing high-virtuality exchanges at any loop order, is 
infrared finite, the soft function collects all non-collinear long-distance singularities, while a 
jet is assigned to each massless parton, collecting the corresponding collinear singularities. 
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This separation is useful because the soft and jet functions which carry the singularities 
have a simple and universal structure, which does not depend on the details of the hard 
process. Most importantly, these functions exponentiate. It is this property that provides 
the basis for resummation, allowing multiple emissions to be derived from a single emission. 

In non-Abelian gauge theories there is a crucial distinction between scattering am- 
plitudes with two coloured partons, and those involving several coloured partons. In the 
former case, only one possible colour flow is present, such that the two partons must nec- 
essarily be coupled by a colour singlet hard interaction. In the latter case, more than 
one possible colour flow may be present. If there are L partons, one may decompose the 
amplitude as 



where the {c^} are a basis of colour tensors, and in the colour index associated with the n**^ 
parton. In this colour flow basis the amplitude is represented by a vector with components 
Aij. Consequently, the factorization between hard and soft modes turns into a matrix 
problem [10, 13-18, 20-35]: the soft function becomes a matrix in colour-flow space, which 
acts on the hard interaction, itself a vector in this space. This adds a degree of combinatoric 
complexity which is absent in the two parton case, and as a consequence the IR singularity 
structure of multiparton amplitudes has been explored only recently relative to that of the 
colour-singlet case. 

In this paper we study soft-gluon exponentiation in multiparton amplitudes. The basic 
tool is the soft (or eikonal) approximation, in which each external hard parton i is replaced 
by a corresponding semi-infinite Wilson line, extending from the origin - where the hard 
interaction takes place - to infinity, in a direction /3j, the direction of motion of parton i. 
The eikonalized amplitude so obtained is much simpler than the original one, and yet it has 
the exact same soft singularities. The simplicity of the eikonal amplitude is due to the fact 
that it only depends on the (Minkowskian) angles between the hard partons and on their 
colours - but neither on their energies, nor their spins. Collinear singularities do depend on 
these, and they are not fully captured by the eikonal approximation. Importantly, however, 
collinear singularities - in contrast to soft ones - are associated with individual massless 
partons: they do not depend on the colour flow in the hard process and can be computed 
using two-parton amplitudes (the Sudakov form factor). Thus, for the purpose of studying 
the exponentiation of soft gluons and addressing the complication posed by colour in the 
multiparton case, we may indeed use the eikonal approximation. In our analysis it will be 
useful to avoid collinear singularities altogether, and we do this by tilting the Wilson lines 
off the lightcone, such that f3f < for any i. At the end one may return to the massless 
case by considering the — >• limit. 

The conventional approach to soft-gluon exponentiation is based on the fact that prod- 
ucts of Wilson lines renormalize multiplicatively [82-85] . Exponentiation is obtained as a 
solution of a renormalization-group equation involving a finite anomalous dimension func- 
tion, the so-called soft anomalous dimension. This approach underlies much of our present 
understanding of the infrared singularity of amplitudes [1-14, 16, 17, 19-38]. Significant 
progress was made in recent years in determining the soft anomalous dimension to two-loop 
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order and furthermore, gaining some insight into its all order properties in the massless 
case [33-35, 37, 38]. It is clear, however, that new techniques are needed to push this 
programme further. 

A complementary approach to exponentiation is the diagrammatic one. This approach 
aims at a direct computation of the exponent of the eikonal amplitude. In the two parton 
case, it has long been established that the exponent may be written as a sum of a restricted 
set of Feynman diagrams - so-called webs - which are irreducible with respect to cutting 
the external parton lines [86-88]^. A direct consequence of irreducibility is that, aside from 
running coupling corrections, each web has only a single pole in dimensional regulariza- 
tion. This goes hand in hand with the fact that the exponent must be expressible as an 
integral over a finite anomalous dimension. Recently, the diagrammatic picture has been 
extended to the multiparton case by two independent groups of authors [89-91]. Both 
groups find that multiparton webs are no longer irreducible, and as a consequence may 
contain higher-order poles, corresponding to subdivergences of the multi-eikonal vertex. It 
is not a priori clear how the singularity structure of multiparton webs can be consistent 
with multiplicative renormalizability and a finite anomalous dimension. 

The main goal of the present paper is to combine the renormalization picture with 
that of webs, in order to gain further insight into the singularity structure of the soft 
gluon exponent in multiparton scattering. We will see in particular that renormalization 
of the eikonal amplitude implies a set of constraint equations for multiparton webs, which 
determine all multiple poles in any given web in terms of lower-order webs. Single poles 
remain the only genuinely new content of webs which is needed to determine the soft 
anomalous dimension at a given order. We will explicitly examine these constraints up to 
four-loop order, and then study how they are realised in a number of non-trivial examples. 
In order to do this, we will introduce a new exponential regulator for disentangling UV and 
IR poles in eikonal calculations at arbitrary orders in perturbation theory. We will also use 
what we learn from this analysis to formulate a new conjecture on the web mixing matrices 
which govern the mixing of colour and kinematic information in multiparton webs [89, 91]. 

The structure of the paper is as follows. In the rest of this introduction, we review 
known results relating to renormalization of products of Wilson lines (section 1.1) and 
multiparton webs (section 1.2), which will be useful in what follows. In section 2 we dis- 
cuss the renormalization of multiparton webs in detail, and formulate constraint equations 
relating the singularities of webs across different orders in perturbation theory, alongside 
expressions for the soft anomalous dimension coefficients in terms of webs. In section 3 
we introduce the infrared regulator we use in order to compute ultraviolet singularities in 
eikonal amplitudes, and then examine how the renormalization constraints of the previous 
section are realised using specific three-loop examples, where certain technical details are 
collected in appendices. We show that for the leading singularity, 0(e~^), the kinematic 
dependence of each diagram factorises into three one- loop kinematic factors, while the 
next-to- leading singularity, 0{e^'^), can be written as a sum of products of lower-order 

^See [40] for a rederivation of these results, which provides a more elegant solution for the exponentiated 
colour factors. 
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diagrams. These properties, in conjunction with the action of the web mixing matrix, pro- 
vide the key to understanding how the renormalization constraints are satisfied. Next, in 
section 4 we show that the factorization of the leading singularity can be generalised to all 
orders for the class of diagrams made exclusively of individual gluon exchanges. Leading 
singularities in different diagrams D in such webs differ only by a combinatorial factor s{D) 
counting the number of orders in which gluon subdiagrams can be sequentially shrunk to 
the origin. In section 5, we formulate a new conjecture regarding web mixing matrices, 
involving a sum rule over column entries, weighted by the above s{D). For the class of 
webs made of single-gluon exchanges this sum-rule directly explains the cancellation of the 
leading singularities required by renormalization. It appears, however, that the sum rule 
itself is completely general, and holds for any web, similarly to other properties of the web 
mixing matrices [89] which already have general combinatorial proofs [91]. In section 6 we 
summarize our conclusions. 

1.1 Wilson lines in multiparton scattering 

As stated in the previous section, the conventional approach to soft-gluon exponentiation 
in scattering amplitudes is based on multiplicative renormalizability of the corresponding 
eikonal amplitude. In this section, we review this subject, preparing for the analysis in 
section 2. Throughout this paper we consider fixed- angle scattering amplitudes with hard 
momenta pi, such that \pi ■ pj\ ^ ^qcD' where all invariants are taken simultaneously 
large. The emission of soft gluons corresponds to the eikonal approximation, in which 
the momenta of the emitted gauge bosons may be neglected with respect to the hard 
momenta pi. In this limit, each external particle may be replaced, as far as the interaction 
with soft gluons is concerned, by the semi-infinite Wilson line 

$ft = Pexp|i5,^"dA/3-A(A/3)|, (1.2) 

where the direction /3j is set by the respective hard parton momenta pi. Similarly, the 
colour representation in which the gauge field is defined is determined by that of parton i. 
The eikonal approximation is useful because it simplifies the calculation while capturing 
all soft singularities. Throughout this paper we consider Wilson lines which are off the 
lightcone, [3f ^ (corresponding to massive external particles), thus avoiding any collinear 
singularities. The soft singularities are then captured by the eikonal amplitude 

5ren. {lij,em) = {0 \^ fS, <^> . ■ . ^ $/3^ | 0)^^^^ (1.3) 

where the colour indices of different Wilson lines remain open. Upon projecting on the 
colour tensors introduced in (1.1), eq. (1.3) becomes a matrix in colour flow space. This 
matrix acts on the hard interaction (finite as e — t- 0) which is a vector in this space: 

M'ipi, em) = 541 iHj, em) ^j(P^) • (1-4) 
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In contrast to the full amplitude or the hard function, which depend on the momenta pj, the 
kinematic dependence of the soft function (1.3) is entirely through the cusp parameters^, 



which, similarly to the eikonal Feynman rules, are invariant under rescaling of the velocities. 
Furthermore the soft function is independent of the spins of the emitting particles. 

Note that the operator in (1.3) has been renormalized: all its ultraviolet singularities, 
including in particular those related to the multi-eikonal vertex, have been removed. Fol- 
lowing Refs. [82-85] we know that it renormalizes multiplicatively, namely there exists a 
factor Z which collects ultraviolet counterterms at each order, and such that 5ren. = S Z 
is ultraviolet finite. Similarly to S, the Z factor is a matrix in colour flow space, so 
their product should be understood as a matrix product (we suppress the colour space 
indices throughout). Clearly, upon renormalization the eikonal amplitude must acquire 
some renormalization-scale dependence. Furthermore, it will be convenient to distinguish 
between the renormalization scale of the coupling /i/j and the renormalization scale n in- 
troduced through the Z factor in the process of renormalizing the multi-eikonal vertex. 

Considered in pure dimensional regularization, all loop corrections in S involve just 
scaleless integrals, and therefore vanish identically. Therefore before renormalization, the 
result is trivial, and one has 



Thus, in dimensional regularization the soft (infrared) singularities we seek to know are 
given by the Z factor. As recognized by different authors [1, 3, 13, 14, 33], this allows one 
to explore soft singularities to all loops by studying the ultraviolet renormalization factor 
of the multi-eikonal vertex in (1.3). 

The analysis of the renormalization of cusp or cross singularities of Wilson lines forms 
the basis of much theoretical work in recent years leading to substantial progress in un- 
derstanding the structure of infrared singularities in multi-leg amplitudes. As we review 
below, these singularities may be encoded in an integral over a so-called 'soft anomalous 
dimension' F. The soft anomalous dimension has been fully determined to two-loop order, 
with both massless [17, 18] and massive partons [23-32, 39]. Moreover, in the massless 
case, stringent all-order constraints were derived [33-35] based on factorization and rescal- 
ing symmetry, leading to a remarkable possibility, namely that all soft singularities in any 
multi-leg amplitude take the form of a sum over colour dipoles formed by any pair of hard 
coloured partons. Despite recent progress [33-38, 92, 93], the basic questions of whether 
the sum-over-dipoles formula receives corrections, and at what loop order, remain so far 
unanswered. Further progress in understanding the singularity structure of multiparton 
scattering amplitudes in both the massless and massive cases requires new techniques to 
facilitate higher-loop computations, or a new viewpoint to give further insight into the 
all-loop structure of the anomalous dimension. 

^The dependence on jij at leading order is often expressed via the corresponding cusp angle, see eq. (3.9). 
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"s(/^|),eiR,/^) = 5uv+iR^ (7ijias(/i|),euv,^) 
= Z {'yij,as{nR),euY,fj) . 



(1.6) 
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Our discussion leading to (1.6) has been quite formal. In order to explicitly compute 
the counterterms entering Z, and as suggested by eq. (1.6), we need to introduce an in- 
frared regulator which disentangles UV and IR singularities, and which does not break the 
symmetries of the problem. With the regulator in place, all poles in e = euv are strictly of 
ultraviolet origin, and using D = 4 — 2e with e > will render the integrals well defined. 
The infrared regulator we use in this paper will be described in section 3.1 below. For our 
general discussion here it is sufficient to assume that such a regularization exists. Denoting 
the energy scale associated with the infrared regulator by m, and the scale associated with 
the renormalization of the multi-eikonal vertex by the regulated version of eq. (1.6) takes 
the form: 

Sren. {jij , Osinn) , n, m) =S (7ij , a, (^|) , €, m) Z {■yij,as{fiji),e,^i) . (1.7) 

We emphasize that while the m-dependent S and 5ren. here are different from their non- 
regularized counterparts in eq. (1.6), the Z factor, which by definition depends on the 
ultraviolet only, is identically the same. By acting on S the Z factor removes its e — )• 
singularities, while making it /i dependent. Although not indicated in (1.7), S^en., while 
finite, usually retains some dependence of e through positive powers. This is precisely the 
situation we will encounter here using a minimal subtraction scheme, where counterterms 
are pure poles. 

To proceed one introduces an anomalous dimension F, defined by 

-ZT. (1.8) 



dlnjj, 



Similarly to Z itself, F is a function of the kinematic variables, 7jj, as well as of as(/i|j) 
and //. Importantly, however, F must be finite, as we now show. In section 2 we will further 
show that in a minimal subtraction scheme, F (unlike 5ren.) is independent of e, aside from 
its dependence through the running coupling; thus the perturbative coefficients F^") at any 
order a" are truly e-independent. 

To show that F must be finite let us first invert (1.7) getting 

5 (e,m) = 5ren. (/^,"i) (e,/i) , (1.9) 

where we suppressed the dependence on jij and asifij^). We note that while both 5ren. 
and Z^^ on the r.h.s. depend on fi, the l.h.s. does not: /i is only introduced through 
the renormalization procedure of the multi-eikonal vertex. Therefore, differentiating with 
respect to fi we have: 

0= Z e,/u) +5ren(/U,m) — -r, • (1.10) 

a m /u am// 



From (1.8) we have: 



dZ dZ-^ „ dZ- 



Z~' = Z =^ =FZ" 

d In dluji d In /i 
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Substituting this into (1-10) and multiplying by Z from the right we obtain: 

= -Sren[fJ',m)T . (1.11) 

am/i 

This implies that T also describes the scale dependence of the renormalized correlator, 
which is by construction free of any ultraviolet singularities. It follows that 

r = -5ren(^, ^ J] 

ain^u 

must itself be ultraviolet finite (finite in 4 dimensions). Since Z is composed of ultraviolet 
poles, the fact that T in (1.8) is ultraviolet finite is a strong constraint on the singularity 
structure of Z. 

The next observation is that the only scales on which T can depend are renormalization 
scales:^ the renormalization scale of the coupling, fiji, and /j, (in particular, given the 
definition (1.8) F cannot depend on the infrared scale m). Therefore, if we choose fiR = /x, 
the only scale dependence that survives in F is through the argument of the coupling. Thus 
we have: 

T = T{ji,,as{^i^e)) , (1.12) 

where we have emphasised that working in = 4 — 2e dimensions, epsilon dependence oc- 
curs through the running coupling = as{pi^, e), which obeys the following renormalization- 
group equation: 

^^^^2 = V + ^oa^ + bial + h2al H ] , ^° ^ 4^ (y '"^^ ~ ' 

^ (1-13) 

Having established the finiteness of the anomalous dimension, the exponentiation of 
soft-gluon singularities is a direct consequence of the defining evolution equation (1.8). 
However, the situation is complicated in multiparton amplitudes by the fact that F is a 
matrix in colour space. This issue was addressed in section 3 of [90] (see also section 
6 in [89]) and will be developed further here, in section 2. Before addressing the general 
multiparton case it is useful to first recall the structure of the exponent in the colour singlet 
case. Equation (1.8) can then be readily integrated to yield an exponential form: 

-^colour singlet = | 2 _/ 2 ^ ^ ('^*-?' ' "'^ ('^^ ' j ' i'^-'^^) 

where we used the initial condition at — )• oo, where the coupling vanishes due to asymp- 
totic freedom. 

Similarly to [94] the singularities in the exponent are generated here through integra- 
tion over the dimensional coupling (note that in our case these are ultraviolet singular- 
ities) . Substituting the expansion of the anomalous dimension 

oo 

F = ^F(")< (1.15) 

n=l 



^Our argument is nothing but the standard 'separation of variables' argument, where the anomalous 
dimension admitting both (1.8) and (1.11) can only depend on the common arguments of the ultraviolet 
object Z and the infrared one 5rcn- 
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it is easy to see that in the conformal case, where bi = 0, namely where dos/dln fi^ = —Og e, 
one obtains: 

^colour singlet|b^=o = | ^ ^ ^^"^ | ' (l-^^) 

Thus, put simply, in the colour singlet case, the exponent (in a conformal theory) has only 
simple poles in e to all orders in perturbation theory. This result, which is straightforward 
to derive, is central to our understanding of the infrared singularity structure of gauge 
theory amplitudes. 

As we show in detail in section 2, upon relaxing the two simplifying assumptions we 
took above, namely including running-coupling corrections, and going to the multiparton 
case, this structure gets modified, giving rise to higher-order poles in the exponent as 
follows: 



explira)„.+ (ir(^)-i?LrW|„ 



2 



+ ih^'"' + 48? [^^'''^'1 - (^or(^) +6,r(i)) + ^^^^^ + } . 

(1.17) 

Indeed, running-coupling corrections generate multiple poles starting from two loops, and 
in addition, due to the non-commuting nature of F^") in the multiparton case, one finds 
multiple pole terms involving commutators starting at three loops. Importantly, however, 
all multiple poles in e, at any order, are determined by lower-order coefficients. This in 
turn implies non-trivial relations between the webs of different orders, relations we shall 
determine in section 2 and explore through examples in section 3. 

1.2 Multiparton vi^ebs 

Having reviewed the use of Wilson lines in multiparton scattering in the previous section, 
in this section we discuss the recently developed diagrammatic approach to soft gluon ex- 
ponentiation, briefly describing the results of [89-91]. In the case of Wilson loops - or 
alternatively amplitudes with a colour-singlet hard interaction - diagrammatic exponen- 
tiation was first studied in the 1980s by Gatheral [86] and by Frenkel and Taylor [87]. 
In this case, as already mentioned briefly above, only irreducible^ diagrams contribute to 
the exponent. By contrast, in the multiparton case, the exponent receives contributions 
from reducible diagrams as well. Refs. [89-91] have shown that these contributions are 
formed through mixing between kinematic factors J~{D) and colour factors C{D') of closed 
sets of diagrams, thus generalizing the concept of webs to the multiparton case. A web 
^(rai,n2,. represents the contribution to the exponent from a set of diagrams which are 
related to each other by permutations of the gluon attachments to each Wilson line, where 
ni is the number of attachments to line /, with I = 1, . . . , L. A multiparton web is then 



*We distinguish between reducible and irreducible (or colour-connected) diagrams. Reducible diagrams 
are ones whose colour factor can be written as a product of the colour factors of its subdiagrams; see e.g. 



91, 95]. 
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given by 

W(n,,n„...,n,) = j;-F(Z))Cp)= ^ ^(D) iiz?D' C(Z)') , (1.18) 

D D,D' 

where C{D) is the so-called Exponentiated Colour Factor (ECF) corresponding to diagram 
D and the sums are over all the diagrams in the set. A general algorithm for computing 
the relevant weh mixing matrices R was developed in [89], and a corresponding closed- 
form combinatoric expression was established in terms of a sum over decompositions of a 
diagram into complete sets of its subdiagrams (eq. (10) in [91]). 

In [89] the web mixing matrices were computed for a range of non-trivial 3 and 4-loop 
examples, revealing interesting mathematical properties: these matrices are idempotent, 
E? = R, and have zero-sum rows, Rdd' = 0, VD. These properties, which were later 
proven to be completely general [91], are crucial for understanding the colour structure of 
the exponent, the interplay between colour and kinematic dependence, and the singularity 
structure. 

The generalization of the familiar properties of colour-singlet webs, namely their "max- 
imal non-Abelian" colour structure and the absence of sub divergences, to the multiparton 
case is non-trivial. Recall that colour-singlet webs are individually irreducible (or colour- 
connected) diagrams. Each diagram that enters the exponent does so with a "maximally 
non-Abelian" ECF. These diagrams have just one overall ultraviolet pole which is asso- 
ciated with the cusp, and any other subdivergence may only contribute to the renormal- 
ization of the strong coupling. Indeed there is a one-to-one correspondence between the 
irreducibility of the colour structure and the absence of subdivergences: as one can easily 
verify graphically, any diagram which has a reducible colour structure can be shown to 
have subdivergences, since its innermost subdiagrams can be sequentially shrunk towards 
the cusp (the hard interaction) without affecting the attachments of the outer gluons. Each 
such shrinking step gives rise to an ultraviolet singularity, a 1/e pole in dimensional reg- 
ularization. The converse is also true, namely any irreducible diagram would have just a 
single overall divergence associated with the cusp. 

On the face of it this picture is entirely broken in the multiparton case. As shown in 
[89-91] only some multiparton webs are individual irreducible diagrams, while most are 
instead formed as in (1.18) through a sum over a set of diagrams, many of which have a 
reducible colour structure, and therefore have subdivergences. One of the most interesting 
observations in [89], however, was that web properties do generalise to the multiparton 
case upon considering the set of diagrams, rather than individual ones. These properties 
are realised through the action of the web mixing matrix. Being idempotent, R acts as 
a projection operator on the space of colour and kinematic factors of the diagrams in the 
set, thus has eigenvalues and 1 (each with nontrivial multiplicity in general). Equa- 
tion (1.18) then states that contributions to the exponents are formed through particular 
linear combinations of kinematic factors multiplying corresponding linear combinations of 
colour factors. The number of such combinations is determined by the rank of R, cor- 
responding to the multiplicity of its unit eigenvalue, while the remaining combinations, 
corresponding to its zero eigenvalue, are projected out. The non-Abelian nature of these 
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[[1,2],[2,1,3],[3]]={3A) [[1,2],[2,3,1],[3]]=(3B) [[1 ,2],[3,2,1 ],[3]]=(3C) 




[[1 ,2],[1 ,2,3],[3]]={3D) [[1 ,2],[1 ,3,2],[3]]={3E) [[1 ,2],[3,1 ,2],[3]]={3F) 



Figure 1. The six 3- loop diagrams forming the 2-3-1 web in which three eikonal hnes are hnked 
by three gluon exchanges. 

colour factors is realised by forming anti-symmetric combinations; in particular, the fully 
symmetric part of the colour factor in each web is projected out owing to the zero-sum 
rows property. This property is expected to hold separately for subsets of diagrams with 
different planarity [91]. 

The pole structure of multiparton webs is rather intricate, but crucially, it is strongly 
constrained by renormalization. The simplest case to analyse is that of the leading singu- 
larities, namely 0(l/e") poles at 0{ag), corresponding to the maximal subdivergence of 
each diagram. As already noted in [89] the only 0(l/e") singularities are the ones purely 
associated with the running coupling (specifically, they are given by eq. (2.16) below). 
Therefore, the maximal subdivergences in webs with two or more connected pieces must 
cancel between different diagrams through the operation of the web mixing matrix. For 
example, in VK(2,3,i) of figure 1 one expects [89] cancellation of 0{e~^) poles among the 
lower three diagrams (diagrams (3A)-(3C) have only subleading divergences), and we shall 
verify these cancellations explicitly in section 3 below. We shall revisit this in section 5, 
where we provide a general explanation of the cancellation mechanism noting a general 
property of the web mixing matrices in the form of a weighted sum over column entries. 

A more complicated structure arises for subleading poles. Starting at ©(l/e""^) at 
0(q"), for n > 3, renormalization of multiparton webs involves commutators of lower- 
orders terms, as in eq. (1.17). This implies that very particular multiple poles should 
survive in multiparton webs. As we shall see below, to be consistent with renormalization, 
these must correspond to commutators between lower-order webs, which are formed from 
subdiagrams of the diagrams in the set. We shall see how this structure arises in section 3. 

In summary, multiparton webs are crucially different to the two parton case in that 
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they are closed sets of diagrams rather than single diagrams. Each individual diagram may 
be reducible, and thus may contain subdivergences associated with the renormalization of 
the multi-eikonal vertex. While in the two parton case all multiple poles in e arise from the 
running of the coupling, in the multiparton case additional multiple poles appear due to 
the non-commutativity of the multi-eikonal vertex counterterms. Nevertheless, similarly 
to the two parton case, all multiple poles at a given order are fully determined by lower 
orders. This highly constrained singularity structure must be matched by diagrammatic 
exponentiation, and therefore by the web mixing matrices. The main aim of the present 
paper is to investigate how this structure is realised. We begin by deriving the singularity 
structure of webs from the renormalization of the multi-eikonal vertex, which is the subject 
of the following section. 

2 Renormalization of multiparton webs 

In the previous section, we reviewed the present state of knowledge regarding soft gluon 
exponentiation in multiparton amplitudes, including the recently developed diagrammatic 
approach. In this section, we will further develop the analysis of the renormalization of the 
multi-eikonal vertex that originated in [90], and was also discussed in [89]. We will see in 
particular that the requirement of multiplicative renormalizability of the eikonal scattering 
amplitude places severe constraints on the singularity structure of multiparton webs. These 
lead to a series of consistency relations expressing the multiple poles of webs at a given 
order in terms of commutators of lower-order web poles. We will develop such equations 
explicitly, and check these in the following section. 

Let us now examine the renormalization of the eikonal vertex along the lines of sec- 
tion 3 in [90] and 6.1 in [89]. As already remarked above, eikonal diagrams are themselves 
scaleless, and thus give zero in dimensional regularization due to the cancellation of UV 
and IR poles. At one-loop order it is relatively straightforward to disentangle these sin- 
gularities. At higher loop orders, the identification of UV and IR contributions becomes 
ambiguous, and in order to compute the UV singularities which comprise the Z factor 
in dimensional regularization, it becomes necessary to implement an additional regulator 
which eliminates all IR divergences. We will indeed introduce such a regulator in section 3, 
where we explicitly set about calculating specific diagrams. With this regulator in place 
all integrals are mathematically well-defined for D < 4 namely for e = (4 — L')/2 > 0. In 
this section, we assume that such a regulator has already been introduced, such that any 
remaining poles are of UV origin. 

As described in the introduction, multiplicative renormalizability of the multi-eikonal 
vertex implies that one may introduce a renormalization factor Z according to eq. (1.7) such 
that Z absorbs all the ultraviolet singularities of 5(e), making the renormalized eikonal 
amplitude 5ren(/f^) finite, and at the same time /^-dependent. Note that both S and Z are 
matrix- valued in colour flow space, consistently with the left-hand side. The ordering of the 
factors on the right-hand side of eq. (1.7) is important because of their non-commuting na- 
ture. This equation constrains the singularity structure that 5(e) may have, thus effectively 
constraining the corresponding webs. 
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Let us see how these constraints arise. Both the eikonal amplitude 5(e) and the Z 
factor exponentiate. That is, one may write [90] 

5(e) = e"', Z(e,/.)=ef, (2.1) 

where the first exponent can be decomposed in perturbation theory as 

oo 

^ = ^^Wa«, (2.2) 

n=\ 

and as (here and below) denotes the renormalized coupling in = 4 — 2e dimensions 
as = as {fJ-'^ , e) . The coefficients on the right-hand side collect all non-renormalized webs at 
a given order in as, i.e. 

w;H < = y 1^/") , , (2.3) 

"1, 

where we have used the notation of section 1.2, adding a superscript to indicate the loop 
order, = 0«). At any given order ©(a") one must have ni + . . . < 2n due 

to the fact that each gluon emission from an external (eikonal) line carries a factor of gs- 
The inequality arises from the fact that gluons may also be connected (by multiple gluon 
vertices or fermion bubbles) away from the external lines, which introduces extra powers 
of the coupling. 

The second exponent in eq. (2.1) may also be decomposed in perturbation theory, as 

oo 

C = J^C(-)a^ (2.4) 

n=l 

where ("^"^ is the counterterm at ©(a"). In general the finite part of the counterterms 
depends on the renormalization scheme. Here we use the MS scheme throughout, meaning 
in particular that C^") does not have any finite parts: it is comprised of only negative 
powers of e. 

As first observed in [90], combining the exponents in eq. (2.1) according to eq. (1.9) 
leads to an exponential form for 5rein where the exponent contains an infinite number of 
commutator terms arising from the Baker-Campbell-Hausdorff formula: 

5ren(At) = exp{ti;(e)} exp{C(e,/i)} 

= exp!^w + C + ^[w,C] + ^{[w, [w, (]] - [C, [w, C]]) - ^ [C, [w, [w, C]]] + •••}. 

(2.5) 

This complication is completely absent in the case of two partons connected by a colour 
singlet interaction, as in that case both 5(e) and Z are not matrix- valued. All commutators 
in eq. (2.5) then vanish, consistent with the fact that, aside from running-coupling terms, 
webs in two-parton scattering have only single poles (a consequence of their irreducibility) . 
In the multiparton case, the commutators are non-zero in general, and lead to a rich 
structure of multiple epsilon poles even in the conformal case. 
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Let us now examine the exponent of eq. (2.5) at the first few orders in perturbation 
theory. In terms of the quantities introduced in eqs. (2.2, 2.4) one finds, up to four loop 
order, 



w 



(1) 



(1) I 



+ ^{2) + 



^(3) = ^(3) + ^(3) + 



W 



(4) 



+ 



W 
1 

12 
1 

24 



(4) + ^(4) ^ 
^(1) 



W 



2 
1 

2 
1 

2 

(1) 



w 



(1)^^(1) 





^(2)" 


+ 




^(1)- 




^(3)" 


+ 




^(2)" 




^(2)- 


+ 







w 



+ 



1 

12 

(3) ^(1) 



(1) _ ^(1) 



W 



+ 



w 



(2) _ ^(2) 



W 



(1) ^(1) 



(2.6a) 
(2.6b) 

(2.6c) 
(2.6d) 



These equations express the sum of renormahzed webs at each order in terms of lower- 
order counterterms and non-renormalized webs. We can already see how this constrains 
multiparton webs: the right-hand side of each equation must be finite as e — )• 0, implying 
the cancellation of all e poles in the particular combinations which occur. The real power of 
these relations, however, is the fact the counterterm is associated with a finite anomalous 
dimension. This implies that C, must have a very particular pole structure, which in turn 
constrains the singularities of the webs. Our next task is therefore to inject this information 
into (2.6). To do this, we begin by defining the anomalous dimension matrix according to 
eq. (1.8) 
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(2.7a) 
(2.7b) 



We have used a known operator identity to relate the derivative of Z = e^ to that of C and 
then the Hadamard lemma to express the result in terms of commutators, see e.g. [96, 97]. 
As for the webs and counterterms, one may expand the anomalous dimension matrix in an 
analogous way, defining 
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n=l 



One may calculate the first few orders explicitly in terms of the expansion coefficients of C 
and the QCD beta function (1.13). To this end we first write 

dos d( 



dC 

dlnfi (iln//2 (^q,^ 

= -2as [e + boas + 6iaf + b2al + 
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(2.9) 



Using (2.7b) we now obtain: 
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(2.10a) 
(2.10b) 
(2.10c) 

(2.10d) 



According to the standard separation of variables argument, hard-soft factorization impHes 
that r, as defined in (2.7a), must be independent of both an infrared cutoff and an ultra- 
violet cutoff, thus it must be finite in four dimensions^ (e — )• 0). Knowing that F^") are 
finite, and working in a minimal subtraction scheme where the counterterms ^^"^ include, 
by definition, only negative powers of e (no ©(e*) for i > 0), we see that r(") are indepen- 
dent of e: neither negative nor positive powers are allowed. It is therefore very convenient 
to express the singularity structure of ("^"^ (and also tt;^"^ below) in terms of r^"'^ 
Inverting the relations in (2.10) we get: 
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(2.11a) 
(2.11b) 
(2.11c) 

(2.11d) 



where we note that owing to a cancellation there are no nested commutator terms in (^^^\ 
Finally substituting the expressions for (^*^") from (2.11) into eqs. (2.6) we extract the 
following results for the singularity structure of the non-renormalized webs w^^^: 
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^Recall that we consider Wilson lines which are off the light cone, so that no coUinear singularities are 
present. 
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In eqs. (2.12) one can see that higher-order powers in e (i.e. not just single poles) are 
indeed present in the non-renormalized webs, as we described in the previous section. 
However, these are fully determined by lower orders. Indeed, as we see below, one can 
explicitly express all multiple poles at any given order in terms of lower-order webs. Finally, 
the new anomalous dimension coefficient F^") will be determined at each order from the 
coefficient of the single 0{l/e) pole in the non-renormalized webs w^^^ of that order. To 
explicitly determine this structure we introduce a further notation for the e expansion of 
(non-renormalized) webs: 
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such that u;^"'"^^ and w^""'^^ are the coefficients of the k^^ negative and positive powers of 
e at C(a"). By inserting (2.13) into (2.12) and extracting the coefficients of for k > 2 
we obtain the relations 
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(2.14a) 
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(2.14c) 
(2.14d) 
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(2.14f) 

These achieve the above objective of encapsulating the constraints derived from renormal- 
ization in a set of consistency conditions involving only coefficients of non-renormalized 
webs. Equations (2.14) can be checked explicitly in specific examples, which forms the 
subject of the following section. 

Before doing this, however, we give explicit expressions for the first few orders of the 
anomalous dimension matrix in terms of the web coefficients of eq. (2.13). These can be 
calculated after extracting the coefficient of in eqs. (2.12), and one finds 

r(i) = -2u;(i'-i) (2.15a) 

r(2) = -4^x;(2'-i) - 2 (2.15b) 

(2.15c) 
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(2.15d) 

As mentioned above, the coefficients F^") themselves contain neither positive nor negative 
powers of e. However, an important observation is that higher-order coefficients in the 
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webs, suppressed by positive powers of e, do in fact enter the expressions for the anomalous 
dimension coefficients. For example, both w^^'^^ and w^^''^^ appear in eq. (2.15). 

Let us now discuss the emerging picture for the pole structure of webs {w) at the 
multi-loop level. According to (2.14a), at two loops only running coupling corrections can 
generate a double pole. In fact, it is true to all-loop orders that the maximal divergence 
y^{n-n) gijjiply given by the one-loop web (single gluon exchange) with running coupling 
insertions: 

ti;^'^'-") = - (-60)"-^ w^^^-^^ , (2.16) 
n 

where 60 is the leading coefficient of the beta function defined in eq. (1.13). In particular, 
as explained in [89], commutator terms cannot contribute to the leading singularity, and 
in a conformal theory we have, at any order. 



= . (2.17) 

fcn=0 



Equations (2.16) and (2.17) can be proven by induction. While we do not present a 
detailed proof, we briefiy explain why it holds, repeating the argument already given in 
section 6.1 of [89]. According to eq. (2.6a) above, the one-loop counterterm is just the 
negative of the simple pole of the one-loop web, 

^(l)=_ly;(l,-l) = lr(l). (2.18) 

e 2e 

Consequently at two loops the commutator term in eq. (2.6b) [u'(i),C^^^] does not con- 
tribute at 0(6-2). Indeed C^^^ in (2.11b) receives no commutator contribution, and the 
0{e~'^) on the r.h.s of (2.6b) must be a pure running coupling term, proportional to ^or'-^^ 
which cancels between ti;^^) and C^^^. Proceeding to three loops, one can easily see using 
eq. (2.6c) how this iterates: none of the commutator terms can contribute to the maximal 
singularity, 0(e^^), and indeed the only 0{e^^) term that appears in C^^) in (2.11c) IS 
proportional to T^^\ with two insertions of 6o- 

With eq. (2.17) at hand it is useful to return to the diagrammatic picture of the webs 
via (2.3), and recall that the contribution of each set of diagrams to if*-"^ takes the form 
of a linear combination of kinematic factors times the corresponding linear combination of 

in) 

colour factors, as determined by the mixing matrix R) 



E <L.,n,)= E E-^(^)<,...,n,)(AI)')^^P'). (2.19) 
n-i,...,nL ni,...,ni D,D' 



Thus, the vanishing of the leading singularity in the conformal limit, eq. (2.17), clearly 
invites the conjecture that 

HD)\oie--) Rt^,,...,njD,D') = 0, yD', (2.20) 

D 

namely that the entries of the mixing matrix are precisely such that the leading poles cancel 
out. This deduction relies upon the fact that the colour factors of individual diagrams 
in the web are linearly independent, and on the assumption that there are no further 
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cancellations between different webs contributing at a given order (n). Although this has 
not been proven in general, one does not expect cancellation of poles between webs. In 
particular, connected multiparton webs readily admit eq. (2.17) on their own, and, as will 
be verified in non-trivial three-loop examples in the next section, so do reducible webs 
composed of single- gluon exchange subdiagrams. 

As we shall see in section 5, owing to the factorization properties of the leading singular- 
ities in a special class of diagrams consisting of individual gluon exchanges, the conjecture 
of (2.20) directly translates into a sum-rule on the columns of the mixing matrix itself, 
which we also verify in a variety of examples. 

Let us now proceed to consider subleading singularities. According to the above anal- 
ysis, at each order n, commutator terms start contributing at ©(e""^). The first such 
contribution is at three loops where a double pole is present in the non-renormalized web 
even when 6„ = 0; in (2.14c) we have: 

u;(3,-2) ^ _1 r (1,-1) (2,-1)1 _ (2.21) 

This simple relation applies to webs that do not include running coupling corrections. We 
will verify it explicitly for two three-loop examples in the following section. 

Considering four-loop webs, one sees directly from eq. (2.12d) that only running cou- 
pling corrections generate poles at 0{e~'^) and ©(e"^), namely 
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While the absence of the leading singularity (2.22a) in the conformal case is expected 
according to (2.17), the absence of the next-to-leading one in (2.22b) looks surprising. 
Indeed, a term such as [ifi,-i, [it)i,-i, it;2,-i]] could have a priori appeared in w^'^~^\ but the 
calculation shows that it is absent. The cancellation of this term appears to be coincidental, 
and does not generalise to higher loops. We have checked explicitly that for five-loop webs^ 
one has 

(2.23) 
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thus showing that 0{e~'^) poles at five loops are present even in the conformal case. 

Finally, the 0{e~'^) terms at four loops, in (2.14f), without running coupling corrections 
read: 
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In this section, we have examined in detail the consequences of the multiplicative 
renormalizability of multi-eikonal vertices, and have derived a number of important re- 
sults. The expressions (2.15) show explicitly how the soft anomalous dimension matrix 



®The expression for the five-loop anomalous dimension coefficient is rather lengthy, however, and we do 
not report this here. 
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may be constructed, up to four loops, from the coefficients of multiparton webs, thus clari- 
fying how the diagrammatic approach to soft gluon exponentiation can be used in practice 
to compute the soft anomalous dimension governing soft singularities in multiparton ampli- 
tudes. Furthermore, we have gained detailed understanding of the singularity structure of 
multipartons webs. We formulated consistency relations which determine all multiple pole 
terms in webs in terms of lower-order webs. We have seen that multiple poles are generated 
both by running coupling contributions and by commutators. The simplest example of the 
latter is given by eq. (2.21) which expresses the double pole of three-parton webs (which 
do not contribute to running coupling corrections) in terms of single poles of lower-order 
webs. In the next section we will examine how this relation is satisfied in the context of 
specific three-loop webs. The ensuing calculations will also allow us to study leading poles 
at arbitrary order for a special class of diagrams in section 4. 

3 The singularity structure of multiparton webs: examples 

In the previous section, we examined in detail the renormalization of multiparton webs, and 
found that it was possible to relate all multiple poles of multiloop webs to poles of lower- 
order webs. The main difference with respect to the colour singlet case is that multiple 
poles remain non-trivial in the conformal limit owing to the presence of commutators. The 
first non-trivial example is the double pole at three loops, where renormalization implies 
the relation of eq. (2.21). 

In this section we focus on two specific examples of three-loop webs, compute their 
singularities, and confirm that these are consistent with the relations of the previous section. 
Beyond illustrating the relations per se, this investigation is important to reaffirm the 
general picture proposed in [89] that multiparton webs, namely the closed sets of diagrams 
whose colour and kinematic factors mix, do indeed share the same properties of connected 
webs. In particular, each such web has on its own, a singularity structure that is consistent 
with renormalization, eq. (2.14) above. Equivalently, no cancellations between webs are 
needed. This is not a priori guaranteed, since webs are defined based on diagrams, and 
as such do not in general form gauge-invariant objects. Discussing the cancellation of the 
leading singularities in (2.19) in the previous section, we have already seen that it is crucial 
to verify that the properties of w^"'\ the sum of all webs at a given order, which follow 
from renormalization, are indeed realised on a web-by-web basis. The question we address 
here is more general: we expect that all the renormalization constraints will be fulfilled by 
individual webs, such that it would be possible to renormalize individual webs. We must 
therefore check that also subleading poles are consistent with the constraints, and to this 
end we will examine in detail how eq. (2.21) is satisfied. 

As we have already mentioned, any calculation of a multiloop eikonal diagram requires 
the use of a regulator for disentangling infrared and ultraviolet singularities. For the 
purposes of this paper, we introduce a method for doing this, which we describe in the 
following subsection. It will likely be useful for future calculations and thus of general 
interest. 
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3.1 Calculating webs: an exponential regulator 

As discussed in eq. (1.6) the calculation of Wilson line amplitudes gives rise to scaleless 
integrals, which vanish in dimensional regularisation. This can be understood as the can- 
cellation of IR and UV singularities, where the coefficients of the latter are relevant for 
determining the renormalization of the multi-eikonal vertex. At one loop, it is possible to 
separate unambiguously the IR and UV poles. At higher loops, this separation can become 
ambiguous. Indeed, previous two-loop calculations of the cusp anomalous dimension [3], 
and 2 — )• 2 eikonal scattering in forward kinematics [10], have employed regulators which 
explicitly suppress IR divergences, such that all remaining poles in dimensional regulari- 
sation are manifestly of UV origin. The advantage of such an approach is clear: one may 
then calculate higher loop diagrams including multiple variable transformations without 
having to keep track of how such transformations mix UV and IR behaviour. Most im- 
portantly, the obtained UV poles, leading as well as subleading, are unambiguous. To this 
end refs. [3, 10] have employed a fictional gluon mass up to two-loop order. 

In this paper we introduce a regulator which is convenient for multiloop, multi-eikonal 
calculations. We work directly in configuration space and introduce an exponential regulator 
along each of the Wilson lines by changing the Feynman rule representing the emission of 
a gluon from a Wilson line as follows: 



fe)/3f 



dX e 



(3.1) 



where the integral is defined in the region where m^J —(if > 0, so that the extra exponen- 
tial factor smoothly cuts off any long distance (A — )• oo) contribution. This prescription 
is guaranteed to yield an infrared-finite result because the rest of the integrand in any 
Feynman diagram can never diverge exponentially for A — )• oo. The exponential damping 
thus removes all IR singularities, such that all remaining poles are of UV (A — )• 0) origin. 
These UV singularities and the corresponding Z factor remain unaffected by the regulator. 

This Feynman rule can be formally obtained from a modified definition of the Wilson 
line of eq. (1.2): 



(m) 

/3» 



P exp <^ i^, 



-m\\ 



(3.2) 



where the exponential damping is now built into the Wilson line which provides the source 
for the soft gluon field. Thus instead of computing the webs with ordinary Wilson lines, we 
compute them starting with a product of exponentially-damped Wilson lines, all involving 
the same infrared scale m: 
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(3.3) 



We emphasise that (3.3) and (3.2), in contrast to (1.3) and (1.2), are not thought of as 
an effective physical description of the interactions of the hard partons with soft gluons. 
Rather the modified Wilson line (3.2) is merely a technical tool that allows one to compute 
the UV singularities that comprise the Z factor one seeks to determine. Feynman diagrams 
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Figure 2. One loop web, where the giuon is emitted between partons i and j, whose kinematic 
part is given by eq. (3.7). 



and webs corresponding to (3.3) will acquire the m dependence through the Feynman rule 
(3.1), but, as already emphasised in (1.7) the Z factor will not. 

The chief advantages of this regulator are twofold. Firstly, there are no problems ap- 
plying it at any loop order. Secondly, the regulator does not break the important symmetry 
of Wilson loop computations, namely invariance under rescaling: 

A ^ kA, ft ^ ft/K . (3.4) 

As an illustration, and in order to introduce notation that will be used later on in 
multiloop diagrams, we consider here the calculation of the one-loop diagram in figure 2 
with this regulator in place. We denote the one-loop diagram with a gluon exchanged 
between legs i and j by 

1,,=T,.T, 4^)(7i,■,/iV"^^e), (3.5) 

where we have explicitly written the colour factor. Using the eikonal Feynman rule of 
eq. (3.1) as well as the configuration-space gluon propagator 

D^,{x) = -Mg,, {-x^y-^ , M = , (3.6) 

the kinematic factor takes the form: 

m?) 2 



F^f (7,, , li^/rr?, e) = ^I'^'gl N ft • /3,- ds ( - (sft - t/3jf^ 

AAr(2e)7ij- [ dxP{x,j^j) 
Jo 

= (5)^ f -^^(26) 7., 2F^ ([1, 1-6], [3/2], ^ + ^ 
where, for later convenience, we defined the propagator-related function. 



(3.7) 



P{x, 7y ) = [x^ + (1 - xf - x{l - x)7ij] ' ^ . (3.8) 

The kinematic dependence is expressed in terms of the cusp parameter 'jij, defined in 
eq. (1.5). A more elaborate example of a three- loop calculation is given in section 3.2.1. 



- 21 - 



In the following sections we shall make use of the leading term in the e expansion of 
eq. (3.7), which we denote by 



9l 
167r2e L 



dx P°(x,7ii) 
2iij coth(^ij; 



(3.9) 



where 



iij = cosh ^(-7ij/2) 



and in the first line P^{x,^ij) is the e — )• limit of P{x,'yij) in (3.8). This is the familiar 
result in the one- loop calculation of the cusp anomalous dimension (see e.g. [3]) where ^ij 
is the Minkowski space cusp angle. 

It is straightforward to expand the hypergeometric function in (3.7) to subleading 
orders in e. Thanks to the IR regulator this expansion is under control, and yields a 
unique result. 

Having described the exponential regulator, we can now set about studying the singu- 
larity structure of specific three-loop webs. This is the subject of the following sections. 



3.2 The three-leg three-loop web W, 



(3) 

(2,3,1) 



In this section we consider the three- loop web composed of the six diagrams in figure 1. 
This web was first studied in [89] , and has a number of interesting properties which make 
it useful as a testing ground for studying renormalization constraints. In particular, not all 
diagrams have a leading divergence: whereas diagrams (3D)-(3F) are 0{e~^), diagrams 
(3^)-(3C) are both 0{e~^), and diagram (3i?) is 0{e~^) i.e. has no sub divergences at all. 
The contribution of this web to the exponent of the eikonal scattering amplitude, including 
the relevant mixing matrix, is [89] 



W, 



(3) 



(2,3,1) 



/^(3^)\ 
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-2 


4\ 
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T{3B) 
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-F(3C) 
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C{3C) 


T{3D) 
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C{3D) 


T{3E) 












-2 
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-2 




C{3E) 


\T{3F)) 
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-2 


1/ 




\C{3F)) 



(3.10) 



We now need to calculate the kinematic parts J-{3A) etc. of each diagram. After combining 
these as in eq. (3.10), we should find that the leading divergence of the web cancels, 
and the remaining terms are such that eq. (2.21) is verified, where the right-hand side of 
that equation contains only those lower-order webs that combine to make the diagrams of 
figure 1. This will then amount to the statement that the renormalization of the exponent is 
consistent on a web-by-web basis, with no mixing between different closed sets of diagrams. 

All of the kinematic integrals we need to compute in this paper can be computed 
using the same general recipe. We thus outline this in detail for one specific example 
in the following subsection and summarise the results for the other kinematic factors in 
appendix A. The details of this method will be useful in section 4. 
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1 2 




Figure 3. Diagram {3D) from figure 1 showing the labels for gluon emissions used in eq. (3.11). 



3.2.1 Calculation of T{3D) 

In order to illustrate how to calculate the kinematic factors needed throughout the paper, 
we consider the explicit example of (31?) from figure 1, which has a leading {0{e~^)) 
divergence. Using the eikonal Feynman rule of eq. (3.1) as well as the configuration space 
gluon propagator (3.6) one finds 



f'OO foo fcyo roo roo foo 

T{2,D) = gli^i^f'N^i^i ■ I32?{h ■ Ps) / dsi / ds2 / dh / dt2 / dt^ / du 

Jo Jo Jo Jo Jo Jo 

^(\/-^(«l+«2) + V-^(*l+*2+t3)+\/-^'^) 



X e 



e-l 



(si/3i - ti/32)'(s2/3i - t2/32)lt3/32 - n/33)"J e[si - S2] e[ti - t2] e[t2 - ts], 

(3.11) 

where /i is the dimensional regularisation scale, and we have labelled the gluon emissions 
according to figure 3. Here the variables Sj, ti and u measure the distance along each 
Wilson line at which gluons are emitted, and we use Heaviside functions to implement the 
ordering constraints implied by the topology of the diagram. Equation (3.11) can be made 
more symmetric by rescaling the distance parameters according to Si — t- Sij \J —0\ etc. to 
give 



T{m = ^(^2)3W37?2723 



X e 



-m(si+S2+tl+<2+t3) 



OO /'OO /*00 /'OO /'OO /'OO 

ds\ \ ds2 / dti / dt2 / dt^ / du 
Jo Jo Jo Jo Jo 

{si + tj- Si ti7i2)(si + tl- S2 i27i2)(ti + - is W723) 



e-l 



X e[si - S2]e[h - t2]e[t2 - ts], 

where we have used again the cusp parameters defined in eq. (1.5). 
As a next step, one may transform according to 



(3.12) 



si 

h 



A 



1 — X I 

after which eq. (3.12) becomes 
J'(3L>) 



S2 
t2 



1- J' 



u 



a 



w 
1 — w 



(3.13) 



^^{^^''f'MH2l2^ I dx [ dy [ dwP{x,ji2)P{y,ji2)P{w,j23) 

o Jo Jo Jo 
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fOO /'oo roo 

X I dX dfi da (Xfia) 
'o Jo Jo 



2e-l ^-m.{X+fi+a) 



X e[Ax - H 0[A(1 -x)- - y)] e[/x(l -y)- aw], 

where P(x,7ij) is defined in eq. (3.8). 

One may proceed further using the additional transformations 



(3.14) 



1 - z 



a 



/3 



1 -r 



(3.15) 



so that eq. (3.14) becomes 



T{3D) = |(^2)3._^3^2^^^^ 



6e-l g-m/3 



X / / dy dw P{x,'-f 12) P{y,-f 12) P{w,-f23) 
Jo Jo Jo 

X f dz[z{l - Z)]2^-1 [ 

Jo Jo 



(3.16) 



X e 



X 



z 



y I- z 



G 



(1 — x) z 
(1 - y) " ~z 



e 



^(i-y) 



1 — r 



The f3 integral gives a singular gamma function r(6e). This represents the overall singular- 
ity obtained by simultaneously shrinking all gluons to the origin in figure 3. The remaining 
integrals over z and r both give additional singularities, as expected from the fact that 
diagram (3D) has sub divergences. One may carry out the r integral by transforming to 
p = r/{l — r), and one finds 



/' 

Jo 



drr 



2e-l 



(1-r) 



4e-l 



G 



41 -y) 



w 

1 (z{l-y) 



1 — r 

2e 



z(l-y)/w 



2e V w 

1 (z{l-y) 
2e 



Fi 2e, 6e; 2e + 1 



z{i-yy 



w 



(3.17) 



w 



2e 



where 2-^1 (o; ^; c; z) is the hypergeometric function. Note that, for the purposes of this 
paper in which we wish to check the relation (2.21) for the first subleading pole, we may 
neglect higher order e terms from the hypergeometric function as shown in the last line of 
eq. (3.17). This leads to a drastic simplification. 

The z integral in eq. (3.16) may be carried out in a similar fashion after transforming 
io C, = z/{l — z), and one ultimately finds 



8e2 



w 



2t 



I f X 1 — X 
mm -, 

V Vy 1 - y 



4e 



P{x, li2)P{y, li2)P{w, 723) (1 + 0{e^)) 



(3.18) 
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where C-\ is the overall normalisation factor: 



„ I p?V" si r(i-e)'r(6e) , , 



In calculating T{3D), eq. (3.18) is already sufficient to express the multiple poles of diagram 
(3-D) in terms of poles of lower-order webs, which is all that is needed to verify eq. (2.21). 

Before proceeding with this analysis, let us summarise the method we have used above, 
which can be generally applied in calculating any multi-eikonal diagram comprising of 
individual gluon exchanges. The general recipe for extracting J-{D) up to and including 
the first subleading pole for a general diagram D belonging to this class is 

1. Insert a factor g-™-*\/-/3^ fQj- ^11 gluon emissions from an external line of 4-velocity /3, 
where t is the appropriate distance parameter. 

2. Represent the ordering of successive emissions along a given line by Heaviside func- 
tions. 

3. Rescale all distance parameters t — )• t/ \/ —P"^- 

4. If si and S2 are parameters associated with the same gluon exchange, transform to 
the variables si = Ax, S2 = A(l — x). Thus, for each separate subdiagram (here a 
single gluon exchange) we have a A-type variable which measures the distance of this 
subdiagram from the hard interaction in units of the infrared regulator 1/m. 

5. Repeat a similar transformation, as in (3.15), for each successive pair of A-type vari- 
ables, until only one remains. The integral over this parameter will give an overall 
gamma function representing the overall UV divergence of the diagram. At each 
stage in this process, one trades two dimensionful parameters representing the dis- 
tance of two separate subdiagrams into one common dimensionful parameter and one 
dimensionless parameter z measuring the relative distance. 

6. Eliminate the Heaviside functions by integrating over the 2:-type variables introduced 
in the previous step. This introduces hypergeometric functions, which can be ex- 
panded in e in order to extract the leading and ffi'st subleading pole of J^{D). 

Using the above prescription, we have calculated the leading and subleading poles of all of 
the diagrams which we need for the rest of the paper. These results are collected in ap- 
pendix A. We shall further apply this technique in section 4 to derive a factorised expression 
for the leading pole of any diagram consisting exclusively of single gluon exchanges. 

This method can also be extended to the next subleading order in epsilon so as to 
compute the contribution of these webs to the soft anomalous dimension coefficient F^^^ 
according to eq. (2.15c). Note that this calculation is technically more involved owing to 
the presence of higher-order terms from the expansion of hypergeometric functions such as 
that of eq. (3.17). 
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(2a) (2b) (2t) (2g) 

Figure 4. The two-loop webs which contribute to H^(2,3,i), labelled as in | 



3.2.2 Expansion of Poles 

In the previous section we showed in detail how to calculate the kinematic parts of all 
web diagrams encountered throughout the paper, collecting other necessary results in ap- 
pendix A. In this section, we expand these expressions, and show that the multiple poles 
of the three-loop web of figure 1 can be expressed in terms of sums of products of poles 
of lower-order webs. These results will then be used in the following section to verify the 
relation (2.21). 

We use a notation for the expansion of the n-loop kinematic factor which is 

analagous to eq. (2.13): 

oo 

^(")(D) = ^ F^^'^^Xd) (3.20) 

k=—n 

Then from eq. (A. 4) one finds 

J-(3'-3)(3^) = ^(3.-3)(3^) ^ F^^'-^\?,C) = 0, (3.21) 

as expected from the fact that these diagrams are not maximally reducible. Furthermore 
one finds 

^(3,-3) (3^) = ^(3,-3) (3^) ^ ^(3,-3) (3^) 

6 2 /■! f3 221 

5.712723 / dxdy du;P°(x,7i2)P°(2/,7i2)^'°(t^,723), 



213 X 37r6e3 







where, as above, P^{x,^ij) is shorthand for P{x,^ij) evaluated at e = 0. 

Comparing this with the leading pole for the one- loop diagram of figure 2, where the 
gluon connects lines i and j, eq. (3.9), we may rewrite eq. (3.22) as 

^(3,-3)(3^) = ^(3,-3)(3^) ^ ^{3.-3)(3i7) = 1 (^^(i'-i)(li2))^ J^^^'^^Hhs)- (3.23) 

This result exemplifies the fact that the leading poles of maximally-reducible diagrams 
can be expressed as products of poles of their irreducible subdiagrams. Furthermore, we 
observe that in this particular case each of the diagrams (3D)-(3-F) has the same numerical 
coefficient in front of the product of one-loop subgraphs. We will see the reason for this in 
section 4. 
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Analogous results are obtained for the two-loop results in this web shown in figure 4 
and labelled as in [89]. The full results are given in appendix A. Their leading pole is at 
0{e^^) and can be written in terms of the poles of one-loop diagrams as follows: 

^(2,-2) (2a) = 7^(2,-2) (25) = 1 t[1'-^^43'^\ (3-24a) 

:F^'^-'\2f) = I (j-fr-^))' , :F('^-'H2g) = 0. (3.24b) 



Note that diagram {2g) has zero contribution at 0{e~'^) as we expect from irreducibility. 
Again we see a factorisation by which higher-order poles can be written in terms of lower- 
order web diagrams, where at leading pole a unique product occurs on the right-hand side. 
By comparing eq. (3.24) with eq. (3.23), one sees that the coefficient of the product of 
lower-order web poles appears to be given by 1/n!, where n is the number of single-gluon 
subdiagrams. We return to this point in section 4. 

Having examined the leading divergence, we now consider the web of figure 1 at the 
first subleading pole ~ 0{e~'^). Looking at the kinematic factors in eq. (A. 4), the next-to- 
leading pole contributions can be divided into three distinct categories, arising from 

1. the expansion of C3; 

2. the expansion of the propagator functions P{x,^ij); 

3. the expansion of the a*^ terms in the integrand, where a is a function of {x, y, z, w}. 

The first two of these are identical for the diagrams in the web, and only contribute at 
0(e~2) for the diagrams which have non-zero poles. They will then be proportional 
to e times that leading pole contribution in each case. The observed cancellation of the 
poles therefore means that the contributions from (1) and (2) above will also cancel when 
the diagrams are combined and we therefore do not consider them further here. 

We instead focus on the contributions from (3) and write these in terms of P^{x,jij) 
and the residue of the pole in C3, 

4-1) = -J^—^. (3.25) 
^ 210 X 37r6 ^ ^ 



The results are: 



^ ' ' 2e2 io V \yJ \^-yJJ (3.26a) 
X e(x - y) p0(x,7i2)P°(2/,7i2)P°(w;,723); 

jp(3,-2)(3^^) = 0; (3.26b) 

,(-1) 7?2723 , , , , „ (x\_ / l-x' 



^(3'-2)(3C) = C(-^) ^ dxdydw (log(^) - log 



X e(x - y) P°(x,7i2)P°(y,7i2)P°(«;,723); 



(3.26c) 
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xp0(x,7i2)P°(y,7i2)P°(«',723); 



X 



(3.26d) 



w JJ (3.26e) 



X p0(x,7i2)P°(y,7i2)P°(u^,723); 

xp0(x,7i2)P°(y,7i2)P°(u^,723). 

We expect the contribution from {3B) at this order to be zero given that the diagram is 
completely irreducible, thus has only a single epsilon pole. As was the case at leading pole, 
one may write the above expressions in terms of poles of lower-order diagrams i.e. in terms 
of the relevant J-^^'^^\lij) and J^^'^~^\2n) contributions. To this end we first consider the 
subleading poles of the two- loop webs, whose leading poles have already been written as 
products of one- loop leading poles in eq. (3.24). Again at two-loops, the contributions at 
next-to-leading pole from the expansion of C2 and P{x,"fij) cancel among themselves, as 
was discussed at three-loops. We therefore concentrate here on the type- (3) contributions 
above, written in terms of P^{x,jij) and the residue of the pole in C2 

4''^ = (3-27) 
The two-loop subleading poles are then 

J-(2'-i)(2a) = 4-'^ lin^l^xdw log(^i-^^pO(x,7i2)P°(«^,723) (3.28a) 

^(2.-i)(26) = 4-^) ^J\xdw log P°(x,7i2)P°(^,723) (3.28b) 

^(2,-i)(2/) = C^-'^ ik j\:,Ay iog(^min(^,i^)) pO(x, 7i2)P°(y, 712) (3.28c) 

(3.28d) 



^(^.-)(2,) = C(-) f f dxd, (log(^) -log([^)) 



X e(x-y)pO(x,7i2)P°(y7i2). 

By inspection, F^'^'^^\2a) = —J^^^'^^\2b) and by explicit calculation, one can also show 
that ^(2,-i)(2/) = -^(2.-i)(25). 

Using these results, we can express the three- loop results at this order in terms of 
diagrams with fewer loops as follows: 

^(3,-2) (3^) = ^-F(2'-i)(25) ^(i'-i)(l23) (3.29a) 
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Figure 5. The two possible decompositions of diagram 3D into a product of a two- loop subdiagram 
times a one-loop one, corresponding to the two terms in (3.29d). Note that the latter equation is 
expressed in terms of {2g), while the actual subdiagram is (2/). However, the two are related by 
J-(2-i)(2/) = -J-(2.-i)(2g). 



jp(3,-2)(3^^ = (3.29b) 

^(3.-2)(3C7) = J-F(2'-i)(25) ^(i'-i)(l23) (3.29c) 

^(3,-2) (3^) = ^ (^^(2,-i)(2a)-F(i'-i)(li2) - 2^(2,-1) (2^)^(1,-1) (i23)^ (3.29d) 

^(3,-2) (3^) = J^(2,-i)(2a)^(i,-i)(i^2) (3.29e) 

^(3,-2)^3^) = ^ (-2J-(2'-i)(2a)-F(i'-i)(li2) - -F(2'-i)(2g)^(i'-i)(l23)) . (3.29f) 

This confirms the statement made above that poles of higher loop webs can be written 
as sums of products of poles of lower-order webs. Unlike the case of the leading pole, 
however, there is in general more than one permissible product on the right-hand side 
of each equation corresponding to the different ways of decomposing each web diagram 
into subdiagrams. This decomposition is illustrated in figure 5 for the case of diagram 
3-D. Furthermore, the coefficients in front of each product of lower-order diagrams are not 
simple factorials, as they were in the leading pole case. Armed with eq. (3.29), we now 
have everything necessary to verify the relation (2.21) for this particular three-loop web. 

3.2.3 Verifying the renormalization constraint 

In the previous section, we used the results collected in appendix A to decompose the 
higher-order poles of the web diagrams in figure 1 in terms of pole coefficients of lower- 
order webs. In this section, we use these results to verify the renormalization constraint 
expressed in eq. (2.21). 

To proceed, one may first expand eq. (3.10) and collect coefficients of each kinematic 
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Figure 6. Illustration of the multiplication rule for colour factors of multiparton diagrams, as given 
in eq. (3.32). 

factor to get 

^(?3,i) = I [3C(3A) - 3C{3C) - 2C(3C) - 2C{3E) + AC{3F)]j^{3A) 

+ ^[- 3C{3A) + 6C{3B) - 30(30) + 0{3D) - 20{3E) + C{3F)\F{3B) 

+ -[- 30{3A) + 30{30) + AO{3D) - 20{3E) - 20{3F)\F{30) 

\ (3.30) 
+ - [0{3D) - 2C(3E) + 0{3F)]F{3D) 

+ \[- 0{3D) + 20{3E) - 0{3F)\F{3E) 

+ ^ [0{3D) - 2C{3E) + C(3F)]^(3F). 

It can be seen immediately after substituting the results of eqs. (3.21, 3.23) that the leading 
poles do indeed cancel to give zero contribution at 0{e~^). That is, this cancellation is 
purely a consequence of the properties of web mixing matrices, and does not depend on 
additional commutators of lower-order webs, which only start to contribute at next-to- 
leading pole. 

We now consider the next-to-leading pole contributions and find by substituting eq. (3.29) 
into eq. (3.30): 

W^/a'-Tn = - l(c{m - 2C(3E) + C(3F))^(2,-i)(2a)^(i,-i)(i^2) 

^"^ y ^ (3.31) 

+ - (C(3A) - 0{30) - 0{3D) + C(3F))^(2,-i)(2^)_^(i,-i)(^^^)_ 

Whilst many cancellations have taken place, a non-zero contribution remains which we 
expect from the renormalisation arguments of section 2. 

We may simplify eq. (3.31) by using the multiplication rule 

C{AB) = 0{A)0{B) (3.32) 

for colour factors of diagrams, where AB is the diagram obtained by ordering subdiagram 
B, then A, outwards from the hard interaction (see e.g. figure 6). This allows us to write 
the colour factors for the three loop web diagrams in terms of the colour factors of lower- 
order diagrams. In particular, for the combinations of colour factors appearing in (3.31), 



-30- 



one has 



C{3D) - 2C{3E) + C{3F) = (C{3D) - C{3E) 



C{li2)[C{2a)-C(2h) 
C(li2), C{2a) - C{2b) 



C{3A) - C{3C) - C{3D) + C{3F) 



C{2g)C{l2z)-C{l2z)C{2g) 

-C(2/)C(l23) + C(l23)C(2/) 

C(l23),C(2/)-C(25) 



C{3E) - C{3F) 

C{2a)-C{2h)]C{li2) 

(3.33a) 
(3.33b) 



We now have the pieces required to independently evaluate the left- and right-hand side of 
eq. (2.21): 

«;(3,-2) = _lL(i,-i), ^(2,-1)1 (3.34) 
f)„=o 6 L J 

for this set of diagrams. We use VF(2,3,i) in place of w when we only consider the contri- 
butions corresponding to this particular subset. The left-hand side comes directly from 
eq. (3.31) 



(3,-2) 
(2,3,1) 



1 

6 L 



C(li2),C(2a)-C(25)l^(2,-i)(2«)^(i,-i)(^^^) 



+ 



C7(l23),C(2/)-C(25)l-F(2.-i)(25)-F(i'-i)(l23). 



(3.35) 



We now turn to the right-hand side. At one- loop, there are two diagrams which will 
contribute to this specific three-loop web: 

= C(li2)-F(l'-l)(li2) + C(l23)-F(l'-')(l23) + ... , 



(3.36) 



where here and in subsequent equations the ellipsis stands for those additional webs which 
do not eventually contribute to the three- loop web of interest. At two- loops, the four 
diagrams in figure 4 are relevant. For each of these the corresponding web is obtained by 
combining the kinematic factor with the corresponding exponentiated colour factors (see 
e.g. eq. (5.65) in [89]) as follows: 

^^(2--i) = i(c(2a) - C{2b)) (^(2,-1) (2a) _ ^(2.-i)(26)' 

+ (C{2g) - C(2/))^(2,-i)(2^) ^ _ _ ^ (3.37) 

'C{2a) - C(26)) J-(2'-i)(2a) + {c{2g) - C{2f)]F^^^-^\2g) + . . . , 



where in the second line we 

haveused J"(2,-i)(2^) = -jr(2,-i)(2a) as follows from eqs. (3.28a) 
and (3.28b). Now substituting these into Eq. (2.21) and restricting to contributions corre- 
sponding to the VF(2,3,i) configuration gives 

.(3,-2) 

(3.38) 



W, 



(2,3,1) 



1 r 

6 
1 

6 



C(li2),(C(2a)-C(26))J^(i'-i)(li2)-F(2.-i)(2a) 
C(l23),(C(25)-C(2/))l^(i'-i)(l23)-F(2'-i)(25). 
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This exactly matches eq. (3.35), and so for this web, the constraints from renormahzation 
have been shown to be satisfied. Note, in particular, that they have been satisfied for this 
web by itself, implying no mixing with other webs. This feature is important in that it 
shows that the notion of individual webs as closed sets of diagrams, first proposed in [89], 
survives after renormahzation of the multi-eikonal vertex. 

In the following section, we illustrate the validity of eq. (2.21) using a second three-loop 



(3) 

(1,2,2,1)- 



example, this time for a web that involves four eikonal lines, W^^ 

(3) 

3.3 The four-leg three-loop web VF^^ ^ 

In the previous sections, we have analysed the web of figure 1 in detail, showing how to 
extract the leading and subleading poles in explicit calculations. We saw that multiple 
poles of web diagrams can be decomposed into sums of products of poles of lower-order 
web diagrams. Furthermore, we explicitly verified that the leading poles, 0{€~^), cancel 
through the action of the web mixing matrix, while the next-to- leading pole, 0{e~^) is 
consistent with the renormahzation constraint of eq. (2.21). 

In this section, we repeat this analysis for a second example, this time considering 

(3) 

a three-loop web involving four eikonal lines, W^(\ 2 2i)' This will provide more evidence 
that the renormahzation constraints of section 2 are indeed realised on a web-by-web basis. 
Given that most of the techniques and results have already been shown in the previous 
section (and in appendix A), we will be much briefer here. 

(3) 

The web W^-^ ^ ^-^ is composed of four diagrams, depicted in figure 7. This web was first 



studied in [89], where its contribution to the exponent of the eikonal scattering amplitude, 
including the relevant web mixing matrix, was shown to be: 



W, 



(3) 

(1,2,2,1) 



/ J-(3a)\ 
T{3b) 
-F(3c) 



T 



( 



\ 



l\ 



/C(3a)\ 
C(36) 
C(3c) 
\C{M)) 



(3.39) 



The kinematic factors for each diagram are given in appendix A, and the leading poles are 




Figure 7. The four 3-loop diagrams forming the 1-2-2-1 web in which four eikonal lines are linked 
by three gluon exchanges, with labels as in [89]. 
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^(3.-3) (3a) = 2^(=^'-3)(36) = 2^(3,-3) (3^) ^ ^(3,-3) (3^) 

The web mixing matrix tells us that these only appear in the combination 

J"(3a) - 2J"(36) - 2J"(3c) + T{3d) 

and therefore the leading poles automatically cancel among themselves, as they must. As 
in the previous case, this means that the contributions at ©(e"^) from pieces common to 
all diagrams, e.g. the r(6e) and the P-functions will also cancel between themselves. We 
therefore do not consider them further here. The remaining contributions to the subleading 
pole are given by 



^(3,-2)(3a) = Cg^M /■'dxdydz P°(x,7i2)P°(y,723)^'°(^,734) 
4e Jo 

X (log(l - x) - log(y) - log(l - y) + log(z)) 
^(3.-2)(35) = Cs^l^xdydz P°(x,7i2)P°(y,723)^°(^,734) 

X (2 log(x) - 2 log(y) + log(l - y) - log(l - z)) 
^(3.-2)(3c) = Cs^l^xdydz P°(x,7i2)^'°(y,723)^°(^,734) 

X (2 log(2;) - 2 log(l - y) + log(y) - log(l - x)) 

J-(3.-2)(3d) = ^3^^ f'dxdydz P°(x,7i2)^'°(y,723)P°(^,734) 
^£ Jo 

X (log(l -y)- log(z) + log(y) - log(l - x)) , 



(3.41a) 



(3.41b) 



(3.41c) 



(3.41d) 



which can again be decomposed in terms of the poles of lower-order diagrams. 

The relevant two-loop diagrams are (2a) and (26) from eq. (3.28) and their mirror 
images about the horizontal. We introduce the new labelling 2ijk where i,j,k label the 
ordering of the legs to which gluons are attached, considered from outside in. So (2a) 
becomes 2i23, (26) is 2321, and their mirror images are respectively 2432 and 2234. The 
correspondence between x, y, z and 'jij can be seen from the arguments of the P-functions. 
This notation introduces some ambiguity since by eqs. (3.28a) and (3.28b) 

-^('•-')(2.,fc) = --F(''-^)(2fe,,) (3.42) 

and so we restrict ourselves to 2i23 and 2432 here. We find 

^{3,-2)(3^) ^ +3 J-(2'-i)(2i23)-^(^'-^)(l34) + 3-^(2.-1) (2432)-F(i'-i)(li2) (3.43a) 

^(3.-2) (35) = +2 ^(2,-i)(2^,3)j-(i,~i)(i3,) - 1^(2.-1) (2432)-F(i'-^)(li2) (3.43b) 

-^(=^■-^^(30) = -l^(2.-l)(2i23)-F(^'-^)(l34) + ^-F(2.-l)(2432)-F(l'-^)(ll2) (3.43c) 

^(3--2)(3d) = _2 ^(2,-i)(2^,3)j-(i,-i)(i3,) _ 2 ^(2.-1) (2432)^(1-1) (I12) . (3.43d) 
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Expanding eq. (3.39) gives 

^£2,2,1) = I (C{3a) - C{3b) - C(3c) + C{3d)^ (-F(3a) - 2T{3b) - 2^(3c) + -F(3(i)) 



(3.44) 



and so 



W, 



(3,-2) 
(1,2,2,1) 



- C(3a) - C{3b) - C(3c) + C{3d) 
6 V 

-F(2'-^)(2l23)-F(l'-l)(l34)+-F(2'-^)(2432)J-(^'-^)(ll2) 



C(l34),C(2i23) - C(232l)l-F(2-l)(2i23)-^^''"'^(l34) 



(3.45) 



C(ll2),C(2432)-C(2234)l-F(2.-l)(2432)-F(l'-l)(li2) 



at this subleading pole order, where the following relationships between the colour factors 
(obtained using eq. (3.32)) have been used: 



C(3a) - C{3b) - C{3c) + C7(3d) = [C(l34), ^(2123) - C(232i) 

= "C(li2),C(2432)-C(2234) 



(3.46) 



Equation (3.45) is the contribution of this web to the left-hand side of equation (2.21). We 
now turn to the right-hand side, 

1 



6 



^(1,-1)^^(2,-1) 



(3.47) 



The relevant one-loop pieces are 

= C7(li2)J-(^'-^)(li2) + C(l34)-F(i'-i)(l34) + ... . (3.48) 

The I23 piece does not contribute because there is no corresponding two-loop web which 
it can combine with^. At two- loops, the diagrams are combined with the following expo- 
nentiated colour factors (see e.g. eq. (5.38) in [89]): 

= ^(C(2l23) - C(232l)) (-F(2'-1)(2i23) " -F^'-^) (232l)) 

+ J(C(2432)-C(2234)) f-F^^'-') (2432) (2234)) + ... 



= (C(2i23) - C(232l))-F(''-^n2l23) + (C(2432) " C(2234))-^(''-'H2432) + . • • , 

(3.49) 

where in the second line we used eq. (3.42). Substituting eqs. (3.48, 3.49) into (3.47) and 
keeping only the parts which correspond to this web then gives 

.(3,-2) _ 1" 



w, 



(1,2,2,1) 



6 L 



C(ll2),C(2432)-C(2234) 



-F(2'-^)(2432)-F(1'-1)(1 



12; 



1 

6 L 



C(l34),C(2i23) - C(232l)l-F(2'-l)(2i23)-^^''"'Hl34), 



(3.50) 



which agrees exactly with equation (3.45). We have thus verified in a second example that 
the renormalization constraints of section 2 are satisfied on a web-by-web basis. 



^There is of course a two-loop diagram where two gluons are exchanged between two distinct pairs of 
legs, [[1], [1], [2], [2]], but this is not a web [89], and thus not part of w^^'~^\ 
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4 Factorization of the leading pole of maximally reducible diagrams 

In the previous section we considered two specific three-loop webs and investigated how 
the renormalization constraints are satisfied. In the webs we considered, each diagram was 
composed exclusively of individual gluon exchanges. We saw that for the leading singu- 
larity of these diagrams, 0{e~^), the kinematic dependence factorises into three one-loop 
kinematic factors, corresponding to sequentially shrinking the three gluon subdiagrams to 
the origin. The next-to-leading sing ularity, ©(e'^), can be written as a sum of products of 
two-loop times one-loop subdiagrams. These factorization properties have proven crucial 
in satisfying the renormalization constraints, namely the cancellation of the leading singu- 
larity and the non-trivial relation (2.21) between three- loop webs and one- and two- loop 
webs. 

In this section we take a step towards an all-order generalization of these results. We 
consider multi-eikonal diagrams at any order but restrict ourselves to ones that 

are composed exclusively of individual gluon exchanges. For this class of diagrams we will 
show that a similar factorization property holds for the leading singularity 0(6""): in a 
given web, for all diagrams that have an (D{e~^) pole, this singularity factorizes into n one- 
loop integrals. Moreover, these leading singularities are all the same up to a combinatorial 
symmetry factor s{D), a non-negative integer number which determines how many different 
ways there are in a given diagram D to sequentially shrink the individual subdiagrams to 
the origin. The quantity s{D) is defined to be zero for those diagrams which do not have a 
leading singularity. A more general definition and some examples for the symmetry factor 
will be given in section 5 (and in appendix B) where we explain its role in the context of 
the web mixing matrix. 

Our first observation is that a web of diagrams which are made of individual gluon 
exchanges satisfies the requirement: 

n 

D k=l 

where the sum is over all the diagrams in the web. This relation, which holds to any order 
in epsilon, may be derived as follows. Each diagram has a kinematic part, after rescaling 
the distance variables associated with each gluon emission according to si — )• \/—f3f etc., 
of the form 

POO fOO f'OO 

/ dSn dti... dtn 

JO Jo Jo 

(s„ + t„ — Sn tnjinjn) 

(4.2) 

Here we have assumed that the fc*^ gluon is exchanged between parton lines ik and j^, and 
is labelled by distance parameters {sk,tk)- The last line denotes the fact that there is a 
product of Heaviside functions, whose role is to order the gluons proceeding outwards from 



(7nji •••7*„i„) / dsi... 



X e 

X 



{sl + t\ - siti7iiji) 
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the hard interaction vertex. This is the only difference between the diagrams {D}, and 
summing over all D gives the complete sum of Heaviside functions 



E[n0](i^^}'i*^»=i' (4.3) 



D 

such that one has 

Y.HD) = ll^f{^,'YMj^,,J dsj dt,e— (^'=+*^)[4 + tl-Sfcifc7.uJ" • (4.4) 
D k=i ^ Jo Jo L J 

This is recognisable as the product of one-loop graphs occurring on the right-hand side of 
eq. (4.1). This result will be useful in what follows. 

Let us now consider a particular maximally reducible diagram D, namely one having 
an ©(e"") pole, thus a non-zero symmetry factor s{D) as defined above. Contributions 
to the leading singularity of this diagram come only from the s{D) different integration 
regions for which the n single-gluon subdiagrams are ordered, and can thus be sequentially 
shrunk to the multi-eikonal vertex (the origin, where the hard interaction takes place). 

To illustrate this point, consider the web diagram (2a) shown in figure 4. We redraw 
this in figure 8, labelling the gluon emissions by Ai, such that the corresponding distance 
parameters for the two ends of the gluon are {si,ti). The Heaviside functions in this 
case impose si < t2, and there are then two distinct integration regions for S2, namely 
■S2 > si and S2 < si. These conditions operate on distance variables on different parton 
lines. However, it is meaningful to compare distances on different lines given that we 



have rescaled distance variables according to si s\l \J etc. The two regions are 
shown schematically in figure 8(a) and (b). Although these diagrams are topologically 
equivalent, they have different properties as the gluon A\ is shrunk to the origin. In the 
first case (s2 > si), A\ may be shrunk to the origin independently of A2, which generates 
a leading singularity (in this case an pole). In the second case, the requirement si < si 
means that as Ai is shrunk to the origin, the end of A2 labelled by S2 is "trapped" at 
the hard interaction vertex, such that the gluon A2 becomes collinear to (52, as shown 
schematically in figure 8(c). This would generate an additional (collinear) singularity were 
the parton lines massless. However, we consider massive parton lines in this paper, so that 
the contribution from figure 8(c) has no further singular behaviour, thus has only a single 
pole, corresponding to the shrinking of Ai. 

Let us now denote by Ai < A2 the integration region in which 

si < S2, t2 , and ti < S2, t2 , (4.5) 

corresponding to the ordered region of figure 8(a). From the above discussion, we see 
that the region Ai < A2 generates the leading pole of the web diagram, but that not 
all subleading poles will be captured, as for the latter, one would also need to include 
configurations such as that shown in figure 8(b). Armed with this notation, we can now 
proceed to analyse the leading pole of a general diagram consisting of single gluon emissions. 

The above argument is straightforwardly generalised to diagrams containing any num- 
ber n of individual gluon exchanges, which contribute at leading pole, 0(e~"). That is, 
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(a) (b) (c) 

Figure 8. The web diagram (2a) of figure 4, where we depict: (a) the integration region S2 > si; 
(b) the integration region S2 < si; (c) the effect of shrinking gluon Ai to the origin when S2 < si- 



the leading pole only receives contributions from the integration region in which all sub- 
diagrams are ordered in the manner of figure 8(a). Denoting the gluon subdiagrams by 
Ak with associated distance parameters Sk and tk along the Wilson lines, a given ordering 
can be written as Aj^^ < < . . . < A.,^^ where vr = (vri,7r2, . . . vTn) is a permutation of 
(1,2,..., n) and where the "<" symbol has the same meaning as described in eq. (4.5). 
We may then write the contribution to diagram D from the ordering vr as 



T{D\tt) 



.k=l 



j j d{sk] d{tk] 



X e 



e-1 



This definition is such that 



(4.6) 



(4.7) 



where the sum is over all orderings vr which are consistent with diagram D. In this way one 
recovers the full leading pole of D by combining ordered regions as defined above. Note, 
however, that the sum does not exactly reproduce the Heaviside functions in eq. (4.2) due 
to the fact that it neglects contributions from integration regions such as that shown in 
figure 8(b). These regions will only contribute at subleading orders in e, hence the 0{e) 
corrections in eq. (4.7). 

We now argue that each such unique ordering gives the same leading pole. First we 
make n paired variable transformations of the form 



such that eq. (4.6) becomes 



F{D\t:) 



m 



2e 




k=l 



Xk 
l-Xk 



1 

ikjk I dxk P{xk,^if.j^.j 





^max 



dXk A? 



(4.8) 



(4.9) 
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where we have scaled Afc — )• Afc/m, and A™^^ is determined from the Heaviside functions, 
depending upon the parameters A^+i, Xk+i and Xk occurring in (4.8). 

The next step is to consider the product of A^ integrals on the right-hand side. These 
may be written as 



n / dXk Xf-^e-^'^ = dXn dXn-i... dXi{XiX2...Xnf'-^e-^^=^^K 

JO Jo Jo Jo 

(4.10) 



k=l 



Here we have used the strict ordering yl^^ < A.,^^ < ■ ■ ■ < ^7r„ , where from dimensional 
considerations each upper limit of integration has been expressed as 

A™*^"" = Afe+i ttfc+i , (4.11) 

where Ok+i is a function of x^+i and Xk, and for later convenience we define ai = 1. 

Equation (4.10) can be computed by making the following change of variables, by 
analogy with (3.15): 

k 

o-fc = y^Afc, Zk-i = —^ (4.12) 

1=1 

where, for convenience, zq = 0. It then follows that the original variables Xk can be 
expressed as 

Afc = O-fc (1 - Zfc_i) = ZkZk+i . . . Zn-l (Tn (1 - Zfc-l) • (4.13) 



The Jacobian is 

0(Al, A2, . . . , A„) ^.,2^3 n-l ^A^A\ 

- Cr2fT3 ...an = Z2Z^Z^ . . . Z^_^ fT„ (4-14) 



d{zi,Z2, . . .,Zn-l,(7n) 

and one also has 

A1A2 . . . A„ = ZiZ^zl . . . zlZ\ (1 - Zi){l - Z2) ... (1 - Zn-l) ^n 

so that eq. (4.10) becomes: 

n / dXk Xf-'e-^'^ = / danal^^-'e-^- 
t_, Jo Jo 



Zk 1 - Zk-1 



(4.15) 



The (Tn integral represents the overall UV divergence corresponding to simultaneously 
shrinking all gluons to the origin. The remaining product of z^ integrals can be sequen- 
tially evaluated starting from Zn-i and working backwards to zi. Each such integral de- 
pends upon a single Heaviside function, and may be carried out using the transformation 
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Pk = — Zk)- Similarly to (3.17) each integration results in a hypergeometric function, 
which becomes trivial upon expansion in e. One ultimately finds 



n i-X^'^^ 2e 4e 2(n-2)t 2(n-l)e 

n / ^1'-=-^' - n^n.) ^ ^ . . . 1^ 1^ (1 + . (4.16) 

Substituting this back into (4.9) yields 
^(Z)|vr) = ( — -2^ 1 -^^^y^|J_7iufcy^ ^^^fc^(^fc'7iuJ («fc(^fc'2;fc_i)) ^ 

(4.17) 

where we explicitly wrote the arguments of at = ak{xk, Xk~i)- At leading pole one may 
ignore subleading contributions from the factors a^^^ and by comparing with (3.7) one 
observes that the leading singularity arising from this ordering is 

HD\n)='-±^llT^^\^.,,,). (4.18) 

k 

Given that there are s{D) such orderings, eq. (4.7) yields the following leading singularity 
of L>: 

^(--")(Z?) = ^n^^-^)(7.,,,). (4.19) 

k=l 

Note that this result is consistent with the expression obtained above for the sum over 
all diagrams D in the set, eq. (4.1), given that the total number of orderings of the n 
subdiagrams is 

^s{D) = nl. (4.20) 

D 

Equation (4.19) is the main result of this section. It shows that the leading pole of a 
maximally reducible diagram composed of single-gluon exchanges can be explicitly related 
to the product of the n constituent one-loop diagrams. The leading poles of different 
diagrams in the web are equal up to the factor s{D) which counts the number of orders of 
sequentially shrinking the different subdiagrams to the origin. In the next section we will 
see that this property, in conjunction with the structure of the web mixing matrix, explains 
the cancellation of the leading poles in these webs, as required by renormalization. 

Before doing so, it is worth noting that there is an interesting constraint on the struc- 
ture of subleading poles, stemming from eq. (4.1). This relation implies that the sum 
over all kinematic parts of a given web gives a complete factorised product of connected 
subdiagrams. Consequently, whereas individual diagrams have subleading poles which de- 
pend upon the ordering of their connected subdiagrams, the sum does not. This in turn 
implies that subleading poles associated with non-connected subdiagrams such as those 
on the r.h.s. of eqs. (3.29) must vanish upon summing over all diagrams in the web. 
One may indeed check that the coefficients of the products J^^'^~^\2A)J^^^'~^\^i2) and 
^■(2.-1) (2(7)7^'^-'^'^-'^^ (723) sum separately to zero after substituting the results into eq. (4.1). 
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Likewise, the coefficients of -F(2,-i)^2i23)J^(^'"^H134) and -F(^'"^H2432)-7^^^'"^H1i2) sum to 
zero in eqs. (3.43). Note that this cancellation is not specific to the particular class of 
diagrams considered here. For webs containing ric connected subdiagrams that may be 
non-trivial (containing three gluon vertices etc.) one may use the appropriate generaliza- 
tion of eq. (4.1), which is 



where is a connected subdiagram with kinematic part J^{Ak)- 
5 A new conjecture for web mixing matrices 

The main aim of this paper has been to investigate the renormalization of multiparton webs 
using the renormalizability of Wilson-line correlators. We have seen how to renormalize 
webs, and learned that a tight singularity structure is realised on a web by web basis. The 
explicit examples considered in section 3 confirm our expectation that the leading pole 
of any given web cancels as a consequence of the properties of web mixing matrices in 
conjunction with the factorization properties of the kinematic factors. We further saw in 
section 4 that the factorization of the leading poles generalises to all orders for a particular 
class of diagrams made of individual gluon exchanges. Since we also expect the cancellation 
of leading singularities to persist to all orders - see eq. (2.20) - it is clear that the cancella- 
tion mechanism we observed at three loops must be general. Building upon this insight, in 
this section we formulate a new conjecture for web mixing matrices which explains these 
cancellations. This conjecture, which we formulate below as a weighted column sum-rule, 
appears to be a general combinatorial property of the web mixing matrices. 

First, we need to precisely define the notion of maximal reducibility. In a given web, 
not all diagrams have leading poles. One example is the web of figure 1, where only three 
of the six diagrams are 0{e~^). A maximally reducible diagram has a number of irreducible 
subdiagrams equal to the number of connected subdiagrams. It is easily checked that this 
condition is verified by diagrams (3-D)-(3F) in figure 1, but not by the first three diagrams: 
(3^) and (3C) have two irreducible subdiagrams and (3-B) has only one. All diagrams, 
though, have three connected pieces, which in this case are all individual gluons. Another 
three-loop example, one with a three-gluon vertex, is that of figure 12. These diagrams have 
just two connected pieces. The leftmost and central diagrams are maximally reducible, as 
the two subdiagrams may be sequentially shrunk to the origin. This is not the case for the 
rightmost diagram. 

Having established the terminology, we may state our conjecture as follows: 

Consider a single web W . For a given diagram D in W , let s{D) be the number of distinct 
ways in which the irreducible subdiagrams of D may be shrunk to the origin, where s{D) = 
if D is not maximally reducible. Then the web mixing matrix satisfies the weighted column 
sum rule 





D 



k=l 





D 
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where the sum is over all the diagrams in W . 



As an illustration of this result, consider the web of figure 1. This web has three 
maximally reducible diagrams: {3D)-(3F). Each of these has three irreducible subdia- 
grams, each of which is a single gluon exchange. Furthermore, there is only one sequence 
in each diagram by which one may shrink all the gluons to the origin: one must start with 
the innermost gluon, and work sequentially outwards. The set of values s{D) for the six 
diagrams may then be written as the row vector 



,(2,3,1) 
'D 



111 



(5.2) 



Equation (5.1) then amounts to the statement 



111 





3 


-3 


-2 


-2 


4\ 




-3 6 


-3 


1 


-2 


1 




-3 


3 


4 


-2 


-2 










1 


-2 


1 










-2 


4 


-2 


v 








1 


-2 


1/ 



0, 



(5.3) 



where we have taken the web mixing matrix from eq. (3.10). A second, and perhaps less 
trivial, example is provided by the web of figure 7, all of whose diagrams are maximally 
reducible. However, whereas diagrams (36) and (3c) have only one ordering by which one 
may shrink the gluons to the origin, in diagrams (3a) and {3d) this ordering is not unique. 
This is because the upper and lower gluons act on completely different parton lines, so that 
the corresponding subdiagrams commute with each other. In both cases one may choose 
either to shrink the upper gluon followed by the lower one, or vice versa. Thus, the vector 
of multiplicity factors s{D) is given in this case by 



Jl,2,2,l) 



2 112 



(5.4) 



and eq. (5.1) for this web is 



2 112 





1 


-1 


-1 


l\ 




-2 


2 


2 


-2 




-2 


2 


2 


-2 


v 


1 


-1 


-1 


1/ 







(5.5) 



Further examples of these symmetry factors may be found in appendix B. 

As justification for this conjecture, let us first consider the special class of diagrams 
considered in section 4, consisting entirely of individual gluon exchanges. There we were 
able to derive a result for the leading pole of any maximally-reducible diagram. We saw 
in eq. (4.19) that the leading pole for any such diagram D reduces to the product of the 
n corresponding one-loop subdiagrams times a symmetry factor s{D). Thus, the leading 
poles of the maximally-reducible diagrams in the set are all equal up to this symmetry 
factor. 
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One may then use the fact, discussed in section 2, that in order to satisfy the renormal- 
ization constraints in the conformal limit, eq. (2.17), leading poles, 0{e~''^), must cancel. 
Assuming that there is no cancellation between different webs (which looks rather unlikely 
given the different kinematic dependence) eq. (2.17) translates into a constraint on indi- 
vidual webs, and given the independence of colour factors of different diagrams in a web, 
this can only be realised if the mixing matrix admits eq. (2.20), 

^ ^("'-")(D)i2j:;_ „^)(D,DO = 0, yD'. (5.6) 

D 

Substituting the leading singularities according to (4.19) into (5.6) then immediately im- 
plies that 

E^(^)<l,...,n,)(AI)') = 0, yD'. (5.7) 

D 

This is the weighted column sum rule of eq. (5.1) which, for this class of diagrams, is 
implied by renormalization. The above argument does not extend, however, to diagrams 
with non-trivial connected subdiagrams involving three or four gluon vertices or fermion 
loops, in which case the number of connected pieces in the diagrams, ric, is smaller than the 
order n. There are two reasons for this. Firstly, the maximal pole obtained by shrinking 
individual subdiagrams to the origin would be in this case 0(e~"'=); therefore this singu- 
larity does not relate to the constraint (2.17) corresponding to ©(e^") in w, but rather 
to renormalization constraints of lower orders, as displayed in (2.14), which necessarily 
involve commutator terms. Secondly, the factorization of the 0{e^^'') pole obtained upon 
shrinking all subdiagrams to the origin for such a diagram takes a more complicated form, 
which does not depend solely on the symmetry factor s{D). 

Nevertheless, it is interesting to examine whether the weighted column sum rule of 
eq. (5.1) still applies, even in the case of webs containing gluon interactions off the eikonal 
lines. Some examples are shown in appendix B, and we find, somewhat remarkably, that 
the sum rule holds in all cases! Given that we presently lack an explanation for this from 
renormalization properties alone, it would be interesting to study the implications of this 
result in more detail. The weighted column sum rule is analagous to other properties of 
web mixing matrices (e.g. idempotence and zero sum rows) whose physical implications are 
better understood, and which have been demonstrated to be general using combinatoric 
methods [89, 91]. Two clear goals for future study are thus to attempt to understand the 
physics embodied by eq. (5.1), and also to seek a combinatoric proof, which may involve 
methods similar to those utilised in [91]. 

Note that even the justification of eq. (5.1) for the special case of diagrams with 
individual gluon emissions falls short of a rigorous proof. The reason for this, as stated 
above, is that we have assumed that no cancellations occur between different webs in writing 
eq. (2.20). This complication would be circumvented upon producing a combinatoric proof 
which relies purely on the properties of individual webs, as has been carried out for the 
zero sum row property in [91]. 

In this section, we have used the insight gained in previous sections from explicit 
calculation of multiparton webs in order to arrive at a general conjecture regarding the 



-42 - 



columns of web mixing matrices. This complements the zero sum row property, and hints 
at the existence of further combinatoric structures underlying web mixing matrices. This 
concludes, however, the analysis which we will undertake in this paper. In the following 
section, we discuss our results and conclude. 

6 Conclusions 

In this paper, we have considered the renormalization properties of multiparton webs, 
taking a further step in developing the diagrammatic approach [89-91] to soft-gluon expo- 
nentiation in multiparton scattering. Webs constitute an efficient way of calculating the 
exponents of eikonal amplitudes directly. They are therefore an important tool in the study 
of multiparton amplitudes, and in particular their infrared singularity structure. 

Building upon [89, 90], we combined in this paper two complementary approaches to 
soft-gluon exponentiation, one which stems from the multiplicative renormalizability of 
Wilson-line correlators and one which utilises a diagrammatic interpretation. 

According to the diagrammatic approach the fundamental objects composing the ex- 
ponent at any given order in perturbation theory are webs. In the multiparton case, these 
are not individual diagrams, but rather closed sets of diagrams that are related to each 
other by permuting the gluon attachments to the Wilson lines. The colour and kinematic 
factors of these diagrams are entangled through mixing matrices. 

The problem of determining soft singularities in amplitudes is equivalent to the prob- 
lem of computing the renormalization of a corresponding eikonal amplitude. Multiplicative 
renormalizability implies that all soft singularities are encodable in a finite anomalous di- 
mension. Owing to recent progress the soft anomalous dimension in multiparton scattering 
is known to two-loop order. Determining its higher-loop corrections and understanding its 
all-order properties is important from both a general field-theory perspective, and a prag- 
matic collider-physics one. The present paper paves the way for such calculations. There 
are several different aspects to this: 

• Firstly, eikonal calculations in dimensional regularization suffer from an inherent 
problem of giving rise to scaleless integrals involving both infrared and ultraviolet 
singularities. Here we provide a prescription to disentangle these singularities at 
any loop order: by using non-lightlike eikonal lines we avoid collinear poles, and by 
introducing an exponential regulator along the Wilson lines (3.1) we eliminate soft 
ones. This facilitates a direct computation of the anomalous dimension associated 
with the multi-eikonal vertex. 

• We have shown that multiplicative renormalizability, and the subsequent finite anoma- 
lous dimension, imply a highly constrained singularity structure for the exponent of 
the eikonal amplitude. Specifically, all multiple poles at any given order are fully 
determined by lower-loop webs. Higher-order poles appear due to two distinct rea- 
sons: running coupling corrections and, in the multi-eikonal case, the non-commuting 
nature of the colour generators. Both these effects are fully understood and are in- 
corporated in the renormalization constraints, which are summarised in eqs. (2.14) 
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through four loops. These relations offer powerful checks on any multiloop eikonal 
computation. 

• Furthermore, we have demonstrated that these renormalization constraints are sat- 
isfied on a web by web basis. Thus, the combinations formed through the operation 
of the web mixing matrices have a singularity structure that is determined by renor- 
malization, where all multiple poles can be computed using lower-order webs. Only 
the single pole terms at each order contain new information, which in turn allows 
one to determine the corresponding contribution to the soft anomalous dimension 
coefficient T^") according to eqs. (2.15). 

• A major difficultly in multiloop computations in dimensional regularization is the 
increasingly high power in e to which the expansion needs to be made. Being able to 
characterize, a priori, the singularity structure of webs, and consequently isolate the 
single-pole contributions to the anomalous dimension facilitates efficient higher-loop 
computations by both analytical and numerical methods. 

We investigated in this paper the way in which the renormalization constraints of 
eqs. (2.14) are satisfied focusing on the class of diagrams composed of individual gluon 
exchanges. This class of diagrams has been convenient in that all divergences in any 
diagram are related to the renormalization of the multi-eikonal vertex, rather than to the 
renormalization of the strong coupling. As a consequence, this analysis is directly valid in 
any gauge theory. 

Our main findings can be summarised as follows: the renormalization constraints are 
realised on a web-by-web basis through the operation of the web mixing matrix, in con- 
junction with the fact that multiple poles in each diagram reduce to sums of products of 
lower-order diagrams. Each term in this sum corresponds to a particular decomposition of 
the original diagram into subdiagrams. At leading pole, any diagram reduces to a product 
of all connected subdiagrams (single gluon exchanges in the case considered here) while 
at subleading poles several terms occur, corresponding to different decompositions of the 
diagram into its subdiagrams, where the latter may be reducible. 

At subleading pole, and starting from three-loop order, the renormalization of the 
multi-eikonal amplitude exponent becomes more complicated owing to the non-commuting 
nature of the colour generators. As a consequence the renormalization constraints include 
commutators of lower-order webs. The first non-trivial relation of this kind, eq. (2.21), 
expresses the next-to-leading pole of three-loop contributions to the exponent in terms of 
a commutator of two-loop and one-loop webs. In section 3 we verified by explicit calcu- 
lations that this relation is satisfied on a web-by web basis. A key ingredient is the fact 
that the next-to-leading poles can indeed be written as a sum over decompositions of any 
three-loop diagrams into two-loop times one-loop subdiagrams. The coefficients in front 
of each product are clearly combinatoric in nature, and are still somewhat mysterious. A 
full analysis of subdivergences would benefit from a more detailed investigation of these 
numbers. 
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The factorization property of the leading singularity has been established in section 4 
for a general n-loop multi-eikonal diagram composed of individual gluon exchanges. We 
found that in a given web, the leading singularities of all maximally reducible diagrams 
are the same up to an overall symmetry factor, s{D), which counts the number of ways 
of shrinking individual subdiagrams to the multi-eikonal vertex. We have seen that the 
cancellation of the leading singularities in this class of webs can be understood through a 
property of the corresponding mixing matrix: the zero column sum rule of eq. (5.1). 

Furthermore, the sum rule (5.1) itself appears to be completely general, working also 
for web diagrams which contain non-trivial connected subdiagrams. It is not yet clear what 
the physics of this result is in terms of renormalization properties. Further insight may be 
gained by corroborating the sum rule via a combinatoric proof, as has been achieved for 
the zero sum row property in [91]. In any case, the weighted column sum rule may act as 
a useful springboard for further investigation of the combinatorics underlying webs. 

In summary, soft-gluon exponentiation can both be derived from the renormalizability 
of Wilson-line correlators and be constructed diagrammatically via webs. Putting these 
two pictures together allowed us to determine the singularity structure of webs and gain 
further insight to the structure of subdivergences. This understanding and the methods we 
have developed pave the way for higher-loop computation of the soft anomalous dimension 
and for a more complete understanding of infrared singularities in scattering amplitudes. 
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A Results for kinematic factors 

In this appendix, we collect results for the kinematic factors of various web diagrams up 
to three loops, which are used throughout the paper in order to verify the renormalization 
constraint of eq. (2.21). All results have been obtained using the method outlined in 
section 3.2.1. 

Firstly, one may consider the one-loop web of figure 2, whose kinematic part is given 
by (3.7), 

Next up are various two-loop webs, shown in figure 4, where we use labels introduced 
in [89]. Their kinematic parts for the leading and first sub-leading poles are given by 

^{2)(2a) = C2^-^ j\xdw (^i-l^yV(x,7i2)P(«),723)(l + 0(e')) (A.2a) 
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j^{2)(^2b) = C2 j\xdw V(x,7i2)P(^,723)(l + 0(e')) (A.2b) 

^(2)(2/) = C2^l\xdy (min(^,[^))'' P(x, 7i2)P(y, 712) (l + 0(e')) 

(A.2c) 

^(2)(2g) = C2^^ j\xdy |^(^^y'-(^l^y'^e(x-y)P(x,7i2)P(y,7i2) 

x(l + 0(e2)), (A.2d) 
where we defined the common factor 

These expressions, and those that follow in eqs. (A. 4) and (A. 6), are valid up to 0{e'^) 
corrections, which is all we require in this study. We note in passing that there is zero 
contribution at 0{e~'^) for diagram (2g) as we expect from the fact that this diagram is 
irreducible. 

In section 3.2 we consider the three-loop web shown in figure 1. The kinematic parts 
of each diagram for leading and first subleading pole are found to be 



X e(x - y) P{x, 7i2)P(y, li2)P{w, 723) (1 + 



(A.4a) 



8e^ ./o 



1 — y \ I I X 



W 



1 — X 

1^ 



1 — f 1 — X 

w J \\yj vi-y 
2 



2e^ 



e(x - y) P{x, 7i2)P(y, li2)P{w, 723) (1 + 0{e^)) 



(A.4b) 



^(3)(3C) = C3%^ f'dxdydw ^^-y 



-4e 



1 - X 



-2e 



-2e^ 



X e(x - y) P(x, 7i2)P(y, li2)P{w, 723) (1 + 



(A.4c) 



^(3) (3^) = C3%^^ tdxdydw ^^-y 



Jo 



w 



2t 



I ( X 1 — X 

mm -, 

V Vy 1 - y 



4e 



X P(x,7i2)P(y,7i2)P(t«,723) (1 + 



(A.4d) 



-46- 



-F(3)(3£;) 



C3 ^^^^"^^^ / dxdydw 



8e2 







I — X 



w 



2e 



I ( X 1 — X 

mm -, 



2e 



1-y 



2e 



mm 



X 1 — X 



X P(x,7i2)P(y,7i2)P(«^,723) (1 + 



4e 



(A.4e) 



X 



(A.4f) 

P{xni2)P{y,li2)P{w,i2-i) (1 + 



where the common factor C3 is as in eq. (3.19): 



In section 3.3 we consider the three-loop web of figure 7. The kinematic parts in this 
instance at leading and first subleading pole are found to be 

r(66) 



-F(3a) = j dxdydz P{x, 7i2)i^(y, l23)P{z, 734) 

/ /I \ 2e /I \2e /I \ 4e\ (A. 6a) 

^(K^) -(^) (^))(-;<^^') ^ 

•^(3^) = C*3^j[|'dxdydzP(x,7i2)P(y, 723)^(^,734) (^) ^ (^) ^ + ^^''^^ 

(A.6b) 

^^'^) = ^^'8f^r'^'^'^^^^'^^^)^^^'^^^)^(^'^ (A.6C) 
X (1 + 0{e')) 

J^{U) = cJ^ j^dxdydzP{x,-ii2)P{y,l2^)P{z,lM) ^ ^ 

X (1 + 0{e')) , 

where 



l-xj (A.6d) 



differs from eq. (3.19) only in the "fij factors. 

B Further examples of symmetry factors 

In section 5 we formulate a weighted column sum rule for web mixing matrices, eq. (5.1). 
This involves the symmetry factors s{D), which for a given diagram D quantify the number 
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[[1], [2, 3,1], [2], [311 [[11, [2,3, 11, [31, [211 [[11, [2,1 ,31, [31, [211 




[[11, [1,2,31, [31, [211 [[11, [2, 1,31, [2], [3]] [[1 J, [1 ,2,3], [21, [311 



Figure 9. Web whose mixing matrix is given by eq. (B.l). 



of ways in which individual connected subdiagrams may be shrunk to the origin. In this 
appendix, we give some further examples of these symmetry factors for various webs to- 
gether with the corresponding mixing matrices. We shall see that in each case the weighted 
column sum rule of eq. (5.1) is admitted. 

Firstly, we consider the web shown in figure 9, whose diagrams are labelled as in [89]. 
The web mixing matrix is 



2 


-1 


-1 


2 


-1 


-1 




[[1], 


[2,3, 1], 


[2], [3]] 


-1 


2 


-1 


-1 


-1 


2 




[[1], 


[2,3,1], 


[3], [2]] 


-1 


-1 


2 


-1 


2 


-1 




[[1], 


[2,1,3], 


[3], [2]] 


2 


-1 


-1 


2 


-1 


-1 




[[1], 


[1,2,3], 


[3], [2]] 


-1 


-1 


2 


-1 


2 


-1 




[[1], 


[2,1,3], 


[2], [3]] 


-1 


2 


-1 


-1 


-1 


2 




.[[1], 


[1,2,3], 


[2], [3]] 



(B.l) 



where we include the vector of labels on the right-hand side so as to define the ordering of 
the diagrams in the matrix. The symmetry factors {s{D)} for this web are 



111111 

It is straightforward to see that with this s{D), the mixing matrix of (B.l) satisfies eq. (5.1). 
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Next, we consider the four loop web of figure 10, whose web mixing matrix is given by 



1 

24 



24 -18 -6 -6 



24 

















-6 -18 -18 
6-6 2 
6 -2 
-2 6 
2 2 
2 
-2 
2 
-2 
-2 
-2 
2 
2 
2 
-2 



-6 -2 
-2 -6 

2 
-2 
-2 
2 
-2 
2 
-2 
2 



18 -18 -6 
-6 -6 -18 
2 2 
2 
6 
2 
6 
2 
2 
2 
2 
2 
2 
2 
2 
2 



-18 -6 12 12 12 12 12 12 
-6 -18 12 12 12 12 12 12 

-2 2 4 4 -4 -4 -4 4 





4 -4 
-4 4 









4 
4 
4 
4 
4 
4 
4 




4 

4 -4 
4 4 



[1,2] 


[3 


1] 


[4 


3], 


[2 


4] 


[1,2] 


[2 


3] 


[3 


4], 


[4 


1] 


[1,2] 


[3 


1] 


[4 


3], 


[4 


2] 


[1,2] 


[2 


3] 


[3 


4], 


[1 


4] 


[1,2] 


[1 


3] 


[3 


4], 


[4 


2] 


[1,2] 


[3 


1] 


[3 


4], 


[2 


4] 


[1,2] 


[3 


2] 


[4 


3], 


[1 


4] 


[1,2] 


[2 


3] 


[4 


3], 


[4 


1] 


[1,2] 


[1 


3] 


[4 


3], 


[2 


4] 


[1,2] 


[3 


2] 


[3 


4], 


[4 


1] 


[1,2] 


[1 


3] 


[3 


4], 


[2 


4] 


[1,2] 


[2 


3] 


[4 


3], 


[1 


4] 


[1,2] 


[3 


1] 


[3 


4], 


[4 


2] 


[1,2] 


[3 


2] 


[4 


3], 


[4 


1] 


[1,2] 


[1 


3] 


[4 


3], 


[4 


2] 


[1,2] 


[3 


2] 


[3 


4], 


[1 


4] 



(B.2) 



The symmetry factors in this case are more comphcated, and are given by 



(2,2.2,2) 



0011111111222244 



Again, one can easily verify that with this s{D), the mixing matrix of (B.2) satisfies 
eq. (5.1). Note that two diagrams are not maximally reducible, as already pointed out 
in [89]. Furthermore, there are 3 topologies of maximally reducible diagram equivalent 
under reflections and rotations, which give rise to 3 possible values for s{D), each with a 
nontrivial multiplicity. 

Next we consider a second four loop example, namely the web of figure 11. In this case 
the mixing matrix is 





6 


-6 


-6 


6 


-6 


6 


6 


-6" 




"[[1], 


[1,2], 


[2,3], [3,4], 


[4]]" 




-6 


6 


6 


-6 


6 


-6 


-6 


6 




[[1], 


[2,1], 


[3, 2], [4, 3] 


[4]] 




-2 


2 


2 


-2 


2 


-2 


-2 


2 




[[1], 


[2,1], 


[2, 3], [3, 4], 


[4]] 


1 


2 


-2 


-2 


2 


-2 


2 


2 


-2 




[[1], 


[2,1], 


[3,2], [3,4] 


[4]] 


6 


-2 


2 


2 


-2 


2 


-2 


-2 


2 




[[1], 


[1,2], 


[2, 3], [4, 3], 


[4]] 




2 


-2 


-2 


2 


-2 


2 


2 


-2 




[[1], 


[1,2], 


[3, 2], [4, 3] 


[4]] 




2 


-2 


-2 


2 


-2 


2 


2 


-2 




[[1], 


[2,1], 


[2, 3], [4, 3], 


[4]] 




-2 


2 


2 


-2 


2 


-2 


-2 


2_ 




.[[1], 


[1,2], 


[3,2], [3,4] 


[4]]. 



(B.3) 
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[[1,2], [3,1], [3,4], [2,4]] 



[[1,2], [2,3], [4,3], [4,1]] 



[[1,2], [3,2], [3,4], [4,1]] 



[[1,2], [2,3], [3,4], [1,4]] 




[[1,2], [3, 2], [4,3], [1,4]] [[1,2], [1,3], [4,3], [4,2]] [[1 ,2], [3,2], [3,4], [1 ,4]] [[1 ,2], [1 ,3], [3,4], [4,2]] 




[[1 ,2], [3,1], [4,3], [4,2]] [[1 ,21,[1 ,3], [4,3], [2,4]] [[1 ,2],[1 ,3], [3, 4], [2,4]] [[1 ,2], [2,3], [4,3], [1 ,4]] 




[[1],[1,2],[2,3],[3,4],[4]1 [[11,[2,1],[3,2],[4,3],[4]] [[1],[2,1],[2,3],[3,4],[4]] [[1], [2,1], [3,2], [3,4], [4]] 




[[1],[1,2],[2,3],[4,3],[4]] [[1],[1,2],[3,2],[4,3],[4]] [[1],[2,1],[2,31,[4,3],[4]] [[1 1,[1 ,21,[3,2],[3,4],[4]] 

Figure 11. Web whose mixing matrix is given by eq. (B.3). 
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and the corresponding symmetry factors are 



(1,2,2,2,1) 
5d 



113 3 3 3 5 5 



There are again three different possible values of s{D). Interestingly in this case, these 
values are all odd numbers. Also here the sum rule of eq. (5.1) is satisfied. 

Let us now consider a few examples of multiparton webs that contain one or more 
three-gluon vertices. A simple one is the three-loop web W^-^ ^ shown in figure 12. The 






[[1],[1,2,2],[2],[]] [[1],[2,2,1],[2],[]] [[1],[2,1,2],[2],[]] 

Figure 12. Diagrams contributing to eq. (B.4). 



mixing matrix is given by 



and the symmetry factors by 



3 


-3 0" 




-3 


3 




_-3 


-3 6. 






(1,3,1) 



[1], [1,2, 2], [2], 
[1], [2, 2,1], [2], 



1 1 



admitting the sum rule of eq. (5.1). 

Next consider the four-loop webs involving two three-gluon vertices, W/ 
figure 13. The mixing matrix in this case is given by: 



(B.4) 



(4) 

(2,2,2)' 



shown in 




[[1,2],[2,1],[1,2],[]] [[1,2],[1,2],[1,2],[1] [[1,2],[2,1],[2,1],[]] 

Figure 13. Diagrams contributing to eq. (B.5). 



[[1,2],[1,2],[2,1],[]] 



-51 - 



1 

24 



and the symmetry factor is 



24 -24 

-24 24 
-24 24 

„(1.3,1) _ 



[[1,2],[2,1],[1,2],[]] 
[[1,2],[1,2],[1,2],[]] 
[[1,2], [2,1], [2, !],[]] 
[[1,2], [1,2], [2, !],[]] 



10 



(B.5) 



Because three out of the four diagrams have no subdivergences, leaving just one maximally 
reducible diagram, [[1,2],[1,2],[1,2],[]], the sum rule of eq. (5.1) reduces to the statement 
that each element in the second row of (5.1) must vanish. This is indeed the case. 
Next, consider the four loop example of figure 14 with mixing matrix 







[[1,21, [1,3,1], [2,3]] 



[[1,2], [2, 3, 2], [1,3]] 



[[1,2], [1,3,1], [3, 2]] 



[[1,2], [2,3, 2], [3,1]] 







[[1,2], [3, 1,1], [3, 2]] 



[[1,2], [3, 2, 2], [3,1]] 



[[1,2], [3, 1,1], [2, 3]] 



[[1,2], [3,2, 2], [1,3]] 







[[1,2], [1,1, 3], [3, 2]] 



[[1,2], [2, 2, 3], [3,1]] 



[[1,2], [1.1, 3], [2, 3]] 



[[1,2], [2,2, 3], [1,3]) 



Figure 14. Diagrams contributing to eq. (B.6). 
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.[[1,2] 


[1,1,3] 


[3,2] 



and symmetry factors 



J2,3,2) 



000000111111 



(B.6) 



(B.7) 



such that the column sum rule is satisfied. Ah of the examples involving three gluon vertices 
that we have so far considered involve s{D) values of 1 or only. As a final non-trivial 
example, it is thus instructive to consider the web of figure 15. Here the mixing matrix is 






[[1], [1,2], [2, 3, 2], [3]] 



[[1], [2,1], [2, 3, 2], [3]] 



aOQQQQQQQQQQC 

[[1], [2,1], [3,2,2], [3]] 






[[1], [1,2], [2, 2, 3], [3]] 



[[1], [2,1], [2, 2, 3], [3]] 



aOQQQQQQQQQQC 

[[1], [1,2], [3, 2, 2], [3]] 



Figure 15. Diagrams contributing to eq. (B.8). 
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_[[1],[1,2],[3,2,2],[3]] 



(Bi 



and the symmetry factors are 



,(1,2,3,1) 



112 2 



(B.9) 



Even in this case, the weighted column sum rule is satisfied. 

Beyond presenting evidence for the sum rule of eq. (5.1), we believe that these ex- 
amples will be useful for further study of the symmetry factors and the role the sum rule 
plays in constraining the corresponding web mixing matrices. It is clear that some inter- 
esting combinatoric patterns underly the possible values of these objects. It would be both 
interesting and useful to know what these are. 
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